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Section 1 


INTRODUCTION 


The purpose of this guide is to train personnel to use the MAESTRO system. The 
MAESTRO system is a mission analysis tool designed to present to the user 
information necessary to make the various decisions required in the design and 
execution of a spaceflight mission. The system was designed so that it can be used 
in both the pre-launch mission planning phase of a mission and during the flight 
as an in-flight decision making tool. The MAESTRO system can be divided into 
the following mission analysis modes : 

1. Trajectory propagation mode 

2. Retro-fire determination mode 

3. Midcourse analysis determination mode 

4. Monte Carlo mode 

5. Verification mode 

6. Orbit Stability mode 

7. Post injection trim mode 

A description of each of the above modes is presented in reference 1. The potential 
user should first read this reference to understand the operation of the mode of 
interest before he attempts use of the program. 

The guide has a section devoted to each program mode. Each of these sections con- 
tains a description of the inputs necessary to run the program mode along with a sample 
case. 

The author hopes that this guide is complete enough so that a new MAESTRO user will 
experience an absolute minimum training period before he is a competent user. Good 
luck with MAESTRO. 
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Section 2 


INPUTS 

Inputs to MAESTRO are made via cards. The input cards are read by MAESTRO's 
input editor which in turn loads the input data into the proper common list. The 
format of the input cards is shown in Table 2.1, "LOC” is the location number of 
the data in the ’’FIELD” portion of the card. Three groups can be placed on each 
card. The value in "LOC” corresponds to the location in the input array where the 
data in "FIELD” is to be stored. A value of LOC greater than 2000 denotes that the 
input corresponds to the gravity model. In this case, the data is not stored in the 
input array. 

A description of the entire input array in shown in Table 2.2. This table describes 
each of the input locations between 1 and 1099 and presents its default value. It is 
envisioned that this table will be used as a reference. Other tables that detail which 
inputs are necessary for the individual modes are presented in the following sections. 

A description of the 2000 series inputs is presented in the next section. 

The input editor has case capability. This capability affords the user a means of 
submitting several runs at one time with a minimum of inputs. Inputs which are 
not changing from run to run should be loaded into the first case. This is the first 
block on input cards. Data changing from run to run are loaded into cases. A 
blank card separates the cases, A run is comprised of all the data read in until 
a blank card is encountered. If the same location number appears more than once, 
the value read in last will be used. For example, if it were desired to fly two trans- 
lunar trajectories which were exactly the same except with different compute intervals, 
then all of the data defining the trajectory and the first compute interval should be 
placed first, followed by a blank card. The second case would follow the blank card. 
This case would only consist of the data defining the second compute interval. A 
submittal can have as many cases as desired. 

There are two important points to note. First, there should not be a blank card at 
the end of the last case. Second, many locations in the input array are initialized 
inside the program. If a pre-initialized location is input via cards, the pre-initialized 
value is lost for all subsequent cases. 
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TABLE 2. 1 


Card Column 
1 

2-6 

7-21 

22-26 

27-41 

42-46 

47-61 


FORMAT OF INPUT CARD 


Description 


Space 

Loe 1, location in the input array, 
right justified 

Data to be entered in the Loc 1 location 
of the input array using D 15.9 format 


Loc 2 ( location in the input array, 
right justified 

Data to be entered in the Loc 2 location 
of the input array using D 15.9 format 

Loc 3, location in the input array, 
right justified 


Data to be entered in the Loc 3 location 
of the input array using D 15 ,9 format 


3 



TABLE 2.2 


MAESTRO INPUT ARRAY* ** 


FORTRAN PRESET 

LOCATION SYMBOL USE*** VALUE DESCRIPTION 


1 


2 

3 

4 


5 


6 

7 

8 
9 


ERRC I 

DELTO I 

DELMN I 

TF I 


Error control for automatic integration 
step size. If = 0 assume in fixed-step 
mode. 

Initial computer step when in automatic 
mode. Note: Need not be input when 
fixed-step. 

Minimum compute interval. Note: Not 
necessary in fixed-step mode. 

Run stop -time. 


PCON I 


Unit conversion factor for positions. The 
value is determined such that, when PCGjf 
multiplies the input units, they will be J 
converted to KMS. The output units will 
be scaled by 1/PCON so that the output 
units will be the same units as input. 


VCON I 1.0 Used the same as PCON except that 

velocities are scaled to KM/SEC. 


EMPTY 

TOL I 

EMPTY 


0.0 


Corrector convergence tolerance in 
twelfth-order predictor-corrector 
integration scheme. Values used depend 
on the accuracy requirements. They 


range from 10"^ to 10 - ^ 


*** 


* unless otherwise stated, the units of the input quantities are: KM, SEC, KG, degrees 

** preset value depends on program operating mode 

I integration V Midcourse verification 

A approach analysis L Lifetime analysis 

M midcourse analysis P PosMnjection trim 
C Monte Carlo analysis 
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LOCATION. 


FORTRAN 
SYMBOL . 


USE 


PRESET 

VALUE 


DESCRIPTION 


10-12 

TMETH(3) 

I 

** 

13-19 

EMPTY 




20 

XMONL 1 




21 

DAYL 




22 

23 

YRL 

HRL 

> 

I 

- 

24 

XMINL 




25 

SECL J 




26-29 

EMPTY 




30 ’ 





31 





32 

33 | 

| ELM(6) 


I 

- 

34 

j 




35 





36 

EMPTY 




37 

DJL 




38 

WTO 


I 

331,40 


39 

ETC 

I 


40 

41 

42 

X(3) 

I 


43 

44 

45 

DX(3) 

I 


46 

DJO 



47 

RA 

I 

- 

48 

DEC 

I 

- 


Times for method table# See KMETH 
(location 1036) for description. 


Launch epoch 


Input initial conditions as orbital 
elements. The order is a,e,f,0J, i, 

The initial conditions may be accepted 
as position and velocity vectors (see 
locations 40-45). The coordinate system 
is defined in KINPT (location 1019). 

Modified julian launch date. Not a program 
input. 

Initial weight 

Ephemeris time correction. If no input, 
the ephemeris time correction will be 
calculated from ETC = 38.66 + .0025921DJ 
where DJ is the number of days since 
the julian date of 2440000.0. 

Initial conditions as position and velocity 
vectors or spherical coordinates. When 
spherical coordinates are used the order 
is velocity, flight path elevation angle, 
flight path azimuth angle, radius, geocentri 
latitude, geocentric longitude. The input 
coordinate system is defined by KINPT 
(location 1019) 

Modified julian launch date. Not a 
program input. 

Initial right ascension and declination 
in Earth equator and equinox of 1950 



LOCATION 


FORTRAN 

SYMBOL 


USE 


PRESET 

VALUE 


DESCRIPTION 


50 

XMONO 

51 

DAY0 

52 

YR0 

53 

HR0 

54 

MIN0 

55 

SEC0 


56-76 

COV 

C 


77-99 

EMPTY 



100-111 

GM(12) 

I 

*** 

112-123 

RE(12) 

I 

*** 

124-135 

WP(12) 

I 

*** 

136-169 

EMPTY 



170-179 

TCOMP(IO) 

I 

** 

180-189 

DELT(10) 

I 

** 


190-192 

XM$#N(3) 

I 

— 

193-195 

DXM$$N(3) 

I 

- 

196 

TMM 

I 

- 

197 

SPRESS 

I 

4.7(10) 


-5 


Initial time 

Input is month, day and year, 
hour, min, second (GMT). 

Upper right-hand triangle of the 
tracking covariance matrix loaded across k 
the rows. Can be input in mean equator g 
and equinox of 1950 or in a local tangent 
plane. Location 1085 determines the system. 
When mean of 1950 is used^the order is | 
X, Y, Z,X, Y, Z. When local tangent plane is 
input ,the order is the position components* 
R, (RxV) x R, R x V and then the velocity | 
components along the same axes. 

Gravitational coefficient. The order of thJ 
planets is defined in locations 1001-1012. 


Planetary radii. 
Planetary rotation rates 


Compute interval table when in fixed comj 
interval mode, {loc (1) = 0.). 

Compute interval =DELT (I) 

when TCOMP (1-1) < T < TCOMP(I) or 

Compute interval = DELT (1) 

when T < TCOMP(l) . Complete compute 

interval table must be loaded when 

alterations are made to preset values. 


Initial position of Moon if osculating 
elements are used for Moon’s position; 


Earth radii and Earth radii/mean solar 



Moon’s epoch in modified julian date. 
Solar pressure at 1 AU, dynes/cm 


*** see Table 3.3 of User's Manual 
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LOCATION 


FORTRAN 

SYMBOL 


USE 


PRESET 

VALUE- 


DESCRIPTION 


198 

199 

REFLK 

EMPTY 

I 

0.2 

Solar pressure reflectivity coefficient 

Tabular thrust history of Motor 1. 
TF1(20) are times since ignition and 

200-219 

TF1(20) 

I, A,M,V 

*** j 

Fl(20) are the thrust values (newtons) 

260-279 

Fl(20) 

I,A,M,V 

*** | 

J 

at the corresponding times. For the 
RAE-B application, this motor is used 
as the midcourse motor. 

220-239 

TF2(20) 

I 

1 

- 1 

j 

> Thrust table for motor 2. 

280-299 

F2(20) 

I 

1 

( 

240-259 

TF3(20) 

I 

1 

- 1 

s Thrust table for motor 3. 

300-319 

F3(20) 

I 

\ 

( 

1 Tabular weight flow rate for motor 1. 

320-329 

TWDOTl(lO) I, A,M, V 

*** 

TWDOT1 is the time past ignition and 
i- WDOT is the flow rate at the corres- 

350-359 

WDOTl(lO) 

I, A,M, V 

*** 

ponding time. The flow rate is in 
KG/sec, 

330-339 

TWDOT2(10) 

I 

- 

Flow rate table for motor 2. 

360-369 

WDOT2(10) 

I 

— 


340-349 

TWDOT3(10) 

I 

- 

Flow rate table for motor 3. 

370-379 

WDOT3(10) 

I 

- 


3 80-382 

TIG (3) 

i,v 

- 

Ignition times for motors 1, 2 and 3. 

383-385 

386-399 

TB0(3) 

EMPTY 

I,V 

- 

Burnout times for motors 1, 2 and 3. 
First and last trial retro firing time on 

400 

TFLREl 

A 

“ 

an approach analysis. Time reference 

401 

TFIRE2 

A 

— 

to liftoff epoch. Used when KAPOPT = 3 

(location 1055) 

402 

DTFIRE 

A 


Increment in retro firing time . 

403 

RAO 

A 

- 

Initial right ascension and declination 
used in the attitude sweep in the approach 

404 

DECO 

A 


analysis. If both zero, velocity vector 
at closest approach is used. 

*** See Table 3. 3 
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FORTRAN 


PRESET 



location 

SYMBOL 

USE 

VALUE 

DESCRIPTION 


405 

DELRA 

A 


Increment in right ascension and 




declination in attitude sweep. 


406 

DELDEC 

A 


Used when KAPOPT = 2 or 3 


407 

EMPTY 





408 

CAR1 

V 

4.0095D8 

Telemetry carrier frequency No. 1 

409 

CAR2 

Y 

4.0D8 

Telemetry carrier frequency No, 2 



A,M, V 

*** 

Observation site geocentric longitude 

410-419 

OBSLON(IO) 

for sites 1-10. 


420 

PSID(l) 

M 

2838. 

Lunar radius 


421 

PSID(2) 

M 

116.5 

inclination w . r . t. target 
planet equator. 

, Midcourse | 

422 

PSID(3) 

M 

. 396D6 

time of flight w. r. t. 

^ analysis 



liftoff 

desired en<J 

423 

PSED(4) 

M 

0.62 

Hyperbolic excess speed 

conditions 

424 

PSID(5) 

M 

0.0 

circular excess speed 


425-429 

PSID(6-10) 

M 


Empty 



DV(3) 



Impulsive velocity of engines 1 to 3. 
Used when KFM0N (location 1047) is set 

430-432 

I 


to 2. Velocities are added at ignition 





times in location (380-382). 






Area of the spacecraft used in solar , 

433 

SOLARA 

I 

13560 

pressure calculation, square centime ter i 

434 

DELTMC 

M 

7200. 

Increment in execution time 


435 

SIGATM 

C,M 

0.7 

One sigma pointing error during 





midcourse maneuver. 


436 

SIGDVM 

C,M 

0.02 

One sigma percentage error of 8 

midcourse velocity. Loaded as percent/1* 

437 

SIGATR 

C 

0.7 

One sigma pointing error during retro. 
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LOCATION 

FORTRAN 
SYMBOL USE 

PRESET 

VALUE 

DESCRIPTION ' 

438 

SIGDVR C 

0.0003 

One sigma percentage error of 
retro velocity. Loaded as percent/100. 

439 

TMC1 C 

36000. 

Time of first correction 

440 

TMC2 C 

324000. 

Time of second correction 

441 

ASPMC A,M,C 

226.0 

Midcourse motor specific impulse 

442 

ASPR A.M.C 

282.5 

Retro Motor specific impulse 

443 

WRETRO A, M, C 

71.44 

Mass of retro fuel 

444 

RO A,M,C,P 

2838. 

Desired lunar orbit radius used in trim 
maneuver 

445 

VC A.M.C.P 


Circular velocity at RO. Program 
calculation, not to be input. 

446 

CIO A, M, C, P 

116.5 

Desired lunar oribt inclination w.r.t. 
lunar equator used in trim maneuver 

447 

448 

BVD(l) M 

BVD(2) M 

6000. 

6000. 

. Midcourse analysis 

Desired B • T desired end conditions. 

Used when IBTR flag 
Desired B • R get tQ ^ (loc 106 2) 

449 

TRINC M,C,P 

0 . 

Tolerance band on inclination correction 
during post injection trim. If inclination 
change is less than TRINC, no inclination 
adjustment is made. 

450-459 

TOUT(IO) I 

TOUT(l) 

1.D20 

= Print table 

460-469 

DTOUT(IO) I 

DTOUT(l)- 

18000. 

Print interval = DTOUT(l) when 
: TOUT(I-l) < T <TOUT(I) or print 
interval = DTOUT(l) when T < TOUT(l) 

470 

ATFULA C 

6.0 

Available attitude control fuel and constant 
to determine attitude control fuel, KG/RAI 

471 

AFUEL C 

0.1678 

Attitude fuel used = AFUEL * attitude 
angle change. 

472 

FT0T C 

20.4 

Total midcourse fuel available 
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FORTRAN 
-location SYMBOL 


USE 


PRESET 

VALUE 


DESCRIPTION 


473 WDROP A,M,C 13.77 


474 CONE A, M 5.0 


475 

TRUE A,M C 

20 

476 

BURNT V 


477 

EMPTY 


478 

TMC M 

7200. 

479 

DINK M 

0.0003 

480-489 

OBSLAT(IO) A,M,V 

*** 


490 

T0L(1) 

M 

10 . 

491 

T0L(2) 

M 

10 . 

492 

T0L(3) 

M 

10 . 

493 

T0L(4) 

M 

.0001 

494 

T0L(5) 

M 

.0001 

495 

T0L(6) 

M 

.02 

496 

T0L(7) 

M 

5. 

497 

T0L(8) 

M 

.2 

498-499 

T0L(9, 10) 



500-511 

RSWTCH(12) 

I 

*** 


Retro drop weight. Weight dropped 
after retro firing. 


Cone angle used in approach analysis. 

The attitude range is this cone angle 
about the velocity vector or input attitude 
when KAPOPT= 1 (loc 1055) . Also used 
in midcourse fixed attitude scan. 

True anomaly range for trial retro firing!# 
in approach analysis, deg. Trial retro 
firings are made from -True to +True 
true anomaly on the approach hyperbola. 
Used when KAPOPT =1. Also used as 
firing range in retro optimization. 

Midcourse motor burn time. Used when 
the intitial state is not obtained from a 
midcourse analysis. 

Initial execution time 
Secant partial step size in midcourse 
analysis also velocity increment when 
using the fixed attitude mode. 

Observation site geocentric latitude 
for sites 1-10 


B* T 
B*R 

time of flight 
hyperbolic excess speed 
circular excess speed 
Total fuel optimization 
closest approach radius 
inclination 

Empty 

Spheres of influence of the planets used 
-to determine the central planet. If the 
distance from the planet is less than theB 
value in RSWTCH the planet is the centrB 
planet. If none of the planets are central* 
then the Sun is considered central. 


Midcourse 
analysis 
tolerances 
desired em 
conditions 


I 

1 


*** See Table 3.3 
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LOCATION 

FORTRAN 

SYMBOL- 

USE 

PRESET 

VALUE 

DESCRIPTION 

512 

RBURN 

A 

20. 

Retro burn time 

513 

TCATST 

I 

0. 

Time to begin closest approach 
testing logic. 

516 

SPNRA 

I 

12.0 

Spin rate, RPM. 

517 

PTI 

I 

238.0 

Initial midcourse motor tank 
pressure, PSIA. 

518 

TPI 

I 

20. 

Initial midcourse motor 
temperature, °C. 

519 

FMC 

I 

20.4 

Current midcourse motor 
fuel, ke. 

520 

FMULT 

I 

1. 

Thrust multiplier. 

521 

WMULT 

I 


Weight multiplier. 

522 

WDOTT 

M 

.0183 

Midcourse motor weight flow used 
with Hamilton Standard thrust in 
midcourse targeting. 

523 

DTM 

M 

0 . 

Compute interva} during motor 
burn in midcourse targeting. 
Compute interval equals burntime 
when value is zero. 
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LOCATION 


FORTRAN 
SYMBOL USE 


PRESET 

VALUE 


DESCRIPTION 


1001 - 

1013 

1014 

1015 

1016 

1017 

1018 


1012 


Bodies in system. Set to 1 if only 
position is necessary, 2 if both 
position and velocity are desired. 
The order is as follows: 


KP(12) I 


1001 MERCURY 

1002 VENUS 

1003 EARTH 

1004 MARS 


1005 JUPITER 

1006 SATURN 

1007 URANUS 

1008 NEPTUNE 


1009 PLUTO 

1010 SUN 

1011 MOON 

1012 ODDBALL 


Note: 1003 =2, 1010 = 1, 
METH 


1011 = 2 are preset 

) Current trajectory propagation method, 
J Not a program input. 


KINT I 


Numerical Integration Scheme 

3. 7th-Order Runge-Kutta 
5. 12th-Order Multistep, 
single step size only 


JL 13 


Launch planet number. 


ENGED A 2 


Engine number of the retro motor. 
Used in a MODE =1 analysis. 


JMN I 5 


KOBLTE I 1 


Lunar and solar ephemeris flag 

1. Mean elements 

2. Mean elements for Sun 
Mean elements + lst-order 
corrections for Moon 

3. Ephemeris tape 

4. Mean elements for Sun, osculating 
elements for Moon (loaded in 
locations 190-196) 

5. Ephemeris tape using Goddard’s 
direct read feature 


Set to 1 for Earth oblateness. Should be 
set to 0 when 1029 = 3. 
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LOCATION 


FORTRAN 

SYMBOL 


PRESET 

VALUE 


USE 


1019 

K3NPT 

I 

1 

1021-1028 

1029 

EMPTY 

KOBL 

I 

* 

1030 

KOUT 

I 

** 

1031 

JT 

I 

11 

1032 

KCRASH 

I 

** 

1033 

1034 

EMPTY 

JOCC 

I 

0 

1035 

MODLEM 

I. L 

1 


1036-1038 KMETH(3) I 
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DESCRIPTION 

Input coordinate system flag 

1. Mean equator and equinox of 1950 

2. Mean equator and equinox of date 

3. Mean ecliptic and equinox of date 

4. True equator and prime meridian 
of launch planet 

5. True equator and equinox of date 

6. True lunar equator and node. 

7. Earth spherical. See location 40-50 
of INPUT common 

Set to the number of the planet in which 
the gravitational field is to be simulated 
using FIELD2. See location 1035. 

If 1, output accoi'ding to print table. If 
zero, output only at beginning and end of 
run. If -1, no output. 

Target planet number 

Closest approach flag 

0 No closest approach test 

1 Continue after closest approach 

2 Stop on closest approach 


Occultation flag. Set to the planet 
number where the observer is located. 

Lunar gravity model flag 

1 for Houston LI model of Lunar field 

2 for Earth J , J , J 

Z O 4 

3 for JPL 15 by 8 Lunar Field 

5 zeroes initial field so new field can be 
input 

10 Used field set in last case 

Trajectory propagation method 

If T < TMETH(1) METH = KMETII(1) 

If TMETH (1-1) <TC <TMETH(I) 

METH « KMETH(I) 

TMETH is in location 10. 

1. Cowell 

2. Encke 

3. NICE/True 

4. NICE/Mean 

5. Averaged Equations of 4 

6. Multiconic, fixed stepsize only 

7. Integrals e cos to, co+ f 

S. Equations 7 are used in averaging 


FORTRAN 

LOCATION SYMBOL 


USE 


PRESET 

VALUE 


DESCRIPTION 


Output coordinate systems - Launch and 
target planets 


1039-1040 

KOUTPT (2) 

I 

3,2 

1. Mean equinox and ecliptic of date 

2. True equator and prime meridian 

3. Mean Earth equator and equinox of 1950 

1041 

JRA 

A, M 

10 

4. True Earth equator and equinox of date 

5. No output 

Number of right ascensions and 
declination angles used in the attitude 
sweep of the approach analysis. Must 

1042 

JDEC 

A, M 

15 

be less than 16. Also used in midcourse 

1043 

KFIRE 

A, M 

20 

fixed attitude analysis. 

Number of trial retro firing times on 

1044 

MODE 

ALL 


approach analysis less than 20 
Program mode 

0 Fly to TF or closest approach 
1, 2 Approach analyses 

3 Midcourse analyses 

4 MontecaiTo analyses 

1045 

KD0P 

V 

0 

5 Midcourse verification mode 

6 Lunar lifetime mode 

7 Post-injection trim 

Set to one for doppler analysis 

1046 

KAPOUT 

A 

0 

Set to one if retro firing time analysis 

is to be output for each attitude 

1047 

KFMOD 

I 

0 

Thrusting mode 
0 Tabular thrust / weight 

1048 

KAPSAD 

A 

1 

2 Impulsive velocity 

3 Hamilton Standard subroutine 

Set to one for lunar orbit shadow 

1049 

KSAD0W 

I 

0 

calculations during approach analysis. 

Shadow calculation flag during trajectory 
propagation: 
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0 No shadows determined. — 

1 Shadow time determined by interpolaticj 
while numerically integrating trajectory. 

2 Osculating orbit used to determine times 
Note: KSADOW=l should not be used when J 


averaging. 



FORTRAN PRESET 

LOCATION SYMBOL USE VALUE DESCRIPTION 

Extra output flag in midcourse analysis. 


1050 MCOUT M 0 


1051 JMC M 10 


1052 KMONTE C 2 


1053 

KMAX 

C 

50 

1054 

KSTART 

C 

17 


1055 KAPOPT A 1 


1056 NGRAPH 


1057 KREAD A,V 


1058 KOUT9 I, L 


1. Outputs targeted solution from PROTO 

2. Prints AV, constraint errors for each 
iteration. 

3. Prints jet iteration information along 
with 2. 

Number of midcourse execution times 
simulated 

Monte Carlo description flag 

1. retro only 

2. midcourse and retro 

3. two midcourses and retro 

If negative, the first correction will be 
calculated from the nominal assuming that 
the tracking data is good. 

Sample size of Monte Carlo analysis. 

Random number starter. 

Approach analysis option flag. 

KAPOPT = 1 Firings made between an 
input range of true anomaly about perigee 
on the approach hyperbola. The attitude 
range is an input cone angle about the 
velocity vector at closest approach. 

KAPOPT = 2 Firings made from asymptote 
to asymptote on the approach hyperbola. 

The attitude range is input as initial right 
ascension and declination, increment in 
right ascension and declination, and number 
of attitude angles. 

KAPOPT = 3 Firings made between input 
times. The attitude range is input as in 2. 

If positive, Element set number, for MAESTR 
graphics data base. 

If, negative, replace (NGRAPH | element set. 

Integer used to obtain initial conditions from 
a midcourse analysis. KREAD corresponds 
to the midcourse correction number desired. 
If KREAD is zero, this option is ignored 
and the initial state must be input. 

Auxiliary peripheral output unit number 
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FORTRAN 
t.OOATIQN SYMBOL 


USE 


PRESET 

VALUE 


DESCRIPTION 


Extra output flag. Used when KOUT9 j^O ' 


1059 

KTERM 

I,L 

0 

. u. Orbital elements printed 

1. Position and velocity vectors 
printed 

2 . Both 0 and 1 

1060 

INPWT 

I 

1 

Input array write flag 

0 no write 

1 write 

1061 

MCUNIT 

A, M, V 

11 

Unit number of auxiliary midcourse 
output unit. 

1062 

IBTR 

M 

2 

Miss vector option flag in Midcourse 
Analysis 

1. Use B*T and B-R loaded in BVD 

2. Use PSID(l) andPSID(2) 





Midcourse guidance law 

1063 

KGLAW 

M 

2 

1. Minimum fuel 

2. Fixed time of arrival 

3. Fixed target energy 

4. Variable target energy 

5. Total fuel optimization 

1064 

NGROPT 

M 

1 

Number of trials for which secant matrix 
is recomputed in Midcourse Analysis. 

1065 

NT 

M 

10 

Number of trials allowed in Midcourse 
convergence. 

1066 

JET 

M 

1 

If set, preliminary targeting will be done 
in the Midcourse Analysis. 

1067 

MCLIM 

M 

100 

Limiting factor in the Midcourse Analysis. 
The midcourse correction is limited to j 


MCLIM * DINK (KM/SEC) on each iterati 
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LOCATION 

FORTRAN 

SYMBOL 

USE 

PRESET 

VALUE 

DESCRIPTION 

1068 

NORD 

I,L 

6 

Number of ordinates in averaging — 
less than 16. 

1069 

INT 

I,L 

3 

Number of intervals in averaging — 
less than 16 

1070 

KPROB 

M 

95 

Output probability for midcourse 
execution error. Second mideourse 
execution time and errors are input 
through locations 440, 435, and 436. 
Negative or zero skips error propagation. 

1071 

IBURN 

M 

6 

Trajectory propagation method during 
finite burn of midcourse motor. Impulsive 
calculations are used when set to .zero. 

1072 

KROUT 

M, C 

0 

Extra output flag during retro optimiza- 
tion calculations. 

1073 

KORECT 

I 

0 

If set nonzero, the derivative at the 
end of an integration step will not be 
calculated. 

1074 

KMTOUT 

C 

0 

Flag used to output initial state for 
each sample in a Monte Carlo analysis. 

1075 

KMETHP 

M 

6 

Trajectory propagation method when 
generating partial derivatives. 

1076 

IFIND 

All 


Element set number of the anchor vector to 
be transferred from the differential 
correction program. 

1077 

KTF 

M 

0 

Flag used to determine the type of 


Midcourse analysis. 

= 0 One-dimensional scan of midcourse 
execution times 

> 0 Two-dimensional scan of midcourse 
execution times and flight times. Flight 
times are scanned in KTF one-hour steps 
beginning at the desired flight time in 
location 422, 

< 0 Two-dimensional scan of midcourse 
execution times and midcourse impulsive 
velocity. The impulsive velocity is 
centered about value, loaded into 426 and 
varied in -KTF steps of size DINK 
(location (479). 
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FORTRAN 
LOCATION SYMBOL 


USE 


PRESET 

VALUE 


DESCRIPTION 


1078 IVTI M 0 


1079 

KHIGH 

I 

0 

1080 

NORMIN 

M, C 

0 


1081 NREV A 0 

1082 KAPWT A 0 


1083 NAPUNT P,L 11 

1084 KSOLP I 0 


1085 KCOV C 0 

1086 

1087 IDATT all 0 


Overburn option key 

0. in-plane retro antiparallel at 

periapsis, 

±1. variable inclination procedure 
approaching above or below desired 
inclination 

± 2, variable periapsis procedure cir- 
cularizing in -plane before or after 
periapsis. 

Farthest approach flag. If set, furthest 
approach will be found. 

Retro optimization flag 

0 in-plane at periapsis maneuver for 
underburns, according to IVTI for 
overburns 

± 1. optimize retro to trim inclination 
in PROTO 

2, Same as 1, but in TARGET also. 

Number of revolutions the transfer orbit 
completes before stopping at closest 
approach in an approach analysis 

Flag used to write information from the 
approach analysis in an auxiliary unit. It 
is set to the output unit number when the 
option is desired. 

Unit numbers where initial conditions 
are stored from a previous approach 
analysis on post- injection trim analysis. 

Solar pressure flag. 

0. no solar pressure 

1. pressure in radial direction only 

(s/c assumed to be sphere) 

2. spacecraft assumed to be a cylinder 
spinning along attitude vector. 

Coordinate system of covariance matrix 

0. mean Equator and equinox of 1950 

1. local tangent plane 

Used in PEST version 

Element set number of attitude file. 

Used when attitude is to be input via 
read from ADP. 
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LOCATION 


FORTRAN 

SYMBOL 


USE 


PRESET 

VALUE 


DESCRIPTION 


10S8 

NATUNT 

all 

12 

Unit number of attitude data. 

1089 

IDSAT 

all 

1234567 

Satellite identification number. 

1090 

INIT 

all 

0 

Initial line number of MAESTRO's 
space allotted on the director's display. 
If INIT is zero, no writes are made on 
the director's display. This option is 
only used during real-time operations. 

1091 

KPIT 

P 

0 

Flag used to indicate post-injection 
trim targeting. 

1092 

ISET 

all 

0 

Element set number when using sub- 
routine PUTELS to transfer state to 
GTDS program. 

1093 

KPLOT 

all 

0 

Set to the plot unit number when 
plotting. 

1094 

NMOD 

I,L 

0 

Maximum number of zonals and 

1095 

NMOD 

I,L 

0 

tesserals used to define the gravity 
field. 

1096 




Used in PEST version. 

1097 

KATMOS 

I 

0 

Atmosphere drag flag. Set to 3 for 
Earth Drag. 

1098 

KASTRT 

all 

0 

Flag used to initiate start-up procedure 
to obtain average elements from 
osculating elements. 

1099 

KTIST 

I 

0 

Flag used to determine thrusters used 
in Hamilton Standard program 





0 both thrusters 

1 first thruster only 

2 second thruster only 
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Section 3 


TRAJECTORY PROPAGATION MODE 

In this mode, MAESTRO is used as a trajectory propagator. Inputs to this mode are 
the most extensive and probably the most important as this mode is also used by 
almost all of the other modes. Table 3. 1 presents the inputs that pertain to this mode. 
Many of the inputs are automatically pre-initialized. The pre-initialized values are 
designed for an Earth-Moon trajectory. Thus, if another type of trajectory is desired, 
manj’’ of the preset values will probably have to be changed. Since the inputs for this 
mode are so extensive, they are divided into catagories according to their use. 

3.1 Initial Conditions 

The initial conditions can be input as position and velocity vectors, orbital elements 
or in an Earth-fixed spherical system. These values are input in locations 30-35 or 
40-45. If the initial semi-major axis in location 30 is zero, the program assumes 
the initial conditions are position and velocity vectors or spherical conditions. If 
location 30 is non-zero, orbital elements are used as the initial conditions even if 
one of the other methods is input. The coordinate system designation is input via the 
flag in location 1019. This flag is also used to indicate the Earth-fixed spherical 
inputs. The time corresponding to the initial state is input in locations 50-55. The 
constants input in locations 5 and 6 are used to convert the units of the initial conditions 
to program units. The program imits are kms and kms/sec. 

3.2 Compute Interval 

The compute interval can be input in one of two ways. The first is via the compute 
interval table in locations 170-189. In this table .constant compute intervals are input 
for segments of the trajectory. The times to switch from one constant compute interval 
to the next are also input. The compute interval can also be determined automatically 
inside the program. When this option is used, the error control limit, initial compute 
interval, and minimum compute interval are input in locations 1 to 3. The error control 
limit determines the accuracy of the propagation. The larger the number.the less 
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accurate. Values of this number between . 001 and . 0001 give reasonably accurate 
results for all the methods except Cowell. When this scheme is used, the error 
control should be between 10 and 100 times larger. There are few restrictions in 
the selection of the compute interval. The twelfth and eighth-order multistep must 
be used with a constant compute interval for the whole flight. Multiconic can also be 
used with the automatic compute interval mode. The use is the same as with numerical 
integration. 

3 . 3 Outputs 

The outputs are obtained with respect to the launch and target planets. The frequency 
of the prints are determined by the print table in locations 450 and 469. In this table, 
constant print intervals are input along with time to switch from one print interval 
to the next. The flags in locations 1030, 1039, 1040 define the frequency and 
coordinate systems of output. Prints can only be obtained at compute points when 
the multiconic propagator is used. 

The trajectory output consists of the time since launch epoch followed by a designation 
of the planet and coordinate system. The third line contains the X, Y, Z components and 
magnitude of the position vector from the designated planet to the spacecraft. The 
next line has the velocity vector and magnitude. The osculating orbit is presented 
next. The order of the elements is semi-major axis, eccentricity, true anomaly, 
argument of perigee, inclination, and longitude of the ascending node. See the 
sample case for an example of the output. 

An alternate abbreviated print is also available. The flags in locations 1058 and 
1059 are used to select this form. The first flag specifies the unit number of 
the prints. Be sure that this input does not conflict with previously defined units. 

The second flag specifies the outputs to be presented-orbital elements, position 

and velocity vectors or both. The flags in locations 1030, 1039, and 1040 also 

apply to this output. When orbital elements are output, the output consists of 

the days since launch epoch, semi-major axis, eccentricity, true anomaly, argument 

of the perigee, inclination^ longitude of the ascending node, the planet number, 

the coordinate system (number corresponds to flag 1039 and 1040), and the periapsis 

radius. When position and velocity vectors are output, the output consists 
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of the hours since launch epoch, position vector, radius planet number, the 
coordinate system (number corresponding to the flag input in 1039 and 1040), 
the velocity vector, and the magnitude of the velocity, 

3.4 Force Model 

A, Planets 

The planets in the system are denoted by the flags in locations 1001-1012. 

A one or two is input in the proper location for the desired planet. If the 
planet will never become the central planet, e.g. , Sun in an Earth -Moon 
trajectory, a one should be input. Otherwise, input a two for the desired 
planets. If the planet is included in the system, the perturbing acceleration 
due to that planet will be determined and included in the force model. 

B. Gravity 

A simplified Earth model containing only J^, J^, will be used by setting 
the flag in 1018 to 1. However, a more extensive gravity model of any 
planet including the Earth can also be used. The flags in locations 1029 
and 1035 are used to select this option. Location 1029 is set to the planet 
number while 1035 is set to the model number. The perturbing acceleration 
will be determined when the central planet is the planet specified in location 
1029, Some gravity models are preset in the program. The flhg in location 
1035 is used to pick up one of these models or to indicate that a new model 
is to be input. 

The following fields are available in MAESTRO: (1) Houston LI model 

of the lunar field, (2) JPL 15-8 model of the lunar field, and (3) Post 

Mariner 1971 model of Mars. A description of these models is shown 

in Table 3,2. The numbers in parentheses represent the corresponding 

model numbers loaded into location 1035, Any changes in these fields 

can be accomplished by additional input through the input array in locations 

starting at 2001, The value for C is loaded into 2000 + j x 16 + (i+1), 

ij 

and the value for S.. is loaded into 2300 + j x 16 + (i+1). For instance, 

C corresponds to 2003, C„ with 2071, S 0 with 2352, etc. 

* • " 6,4 3,3 
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The C and S arrays are the cosine and sine terms, respectively, of the 
gravity field expansion. For a further description of these arrays, see 
reference 1, The elements input in the 2000 series override elements 
preset in the chosen model. Thus, to change one element in the JPL 15-8 
field, MODLEM (1035) should be set to 3 and the new element should be 
input in the proper 2000 location. When new fields are to be input, 

MODLEM should be set to 5 to zero out the initial field. If the field used in 
the previous case is to be used again with minor changes, then MODLEM 
should be set to 10 and the appropriate changes input in the 2000 section. 

If a new field is to be input, or the size of a preset model is to be changed 
through inputs, then the flags INMOD{1094) and IMMOD(1095) should be 
set to the new size. INMOD represents the number of zonals and IMMOD 
represents the size of the tesseral array. For instance, if the JPL 15-8 
field is to be expanded to a 16-9 field by adding a zonal and several 
te s serai s , the INMOD should be set to 16 and IMMOD should be set to 9, 

It is important to note that the 2000 section inputs must be the last inputs 
in each case and that these inputs must always be on separate cards 
from the 1 to 1099 inputs. 

The acceleration due to the field is only determined when the planet 
number input in location 1029 is the current central planet. Care should 
be taken not to use the Earth field twice. If an Earth model is used and 
1029 set to 3, be sure location 1018 is zero. 

C, Thrusting 

The program has the capability of firing three separate engines; however, 
only one engine can be burning at a time. The engine thrust characteristics 
are loaded into the tables in locations 200 to 319. In these tables ,the thrust 
values are input at the corresponding times since ignition. Linear interpolation 
is used to obtain the thrust at times between the input values. The weight 
flow rate is input in locations 320 to 379. In these tables, the weight flow 
rate is input at the corresponding times since ignition. Linear interpolation 
is used to obtain the flow rate at times between the input times. The 
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program integrates the weight flow rate to obtain the current spacecraft 
weight. The ignition and burnout times of the engines are input in location 
380-385. The engine may also be simulated by an impulsive velocity 
approximation. When this approximation is used, the impulsive velocity 
is input in location 430 to 432 and the flag in 1047 must be set. The 
ignition times in locations 380-382 are used to determine when the velocities 
are applied. 

D. Solar Pressure 

The acceleration due to solar pressure is included when the flag in location 
1084 is non-zero. When this option is used, the solar pressure at one 
astronomical unit and the spacecraft area and reflectivity are also required. 
These constants are input in locations 197, 433, and 198, respectively. 

The solar pressure force can be determined in two different ways. One 
of the methods considers the spacecraft to be a sphere. Thus, the solar 
pressure force is only acting along the radial rays from the Sun. The 
other method assumes that the spacecraft is a spinning cylinder with the 
spin axis along the spacecraft attitude, input in locations 47 and 48. In 
this model, the solar pressure force consists of two components. One, 

Fr, in the radial direction from the Sun and the other, Fn, normal to the 
spin axis, see the figure below. 



The flag in location 1084 is used to select the desired solar pressure 
model. 
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3.5 Phys Leal C onstants 

Although all the physical constants are pre-set in the program, the user has the 
option to change any of these values. A table describing the preset constants 
is shown in Table 3.3. 


3.6 Program Options 

There are a variety of program options keyed to input flags. The table below 
describes the options available to the user. 


Input Location 
*1014 

*1017 

1018, 1029, 1035 


*1019 

1032 

*1036-1038 


Description 

This flag is used to select the numerical integration scheme 
desired when the state is to be propagated by numerical 
integration. The 12th-order multistep method can only be 
used with a constant compute interval. 

This flag is used to select the desired ephemeris. 

These flags are used to describe the desired gravity 
mode. The 1018 flag denotes a simplified Earth model. 

If a more complex model for the Earth or any planet is 
desired, the flags in locations 1029 and 1035 are to be 
used. See Table 3.2 for a description of the available 
preset gravitational models. 

This flag is used to indicate the coordinate system of 
the input spacecraft state. 

This flag is used to indicate the stopping criterion of 
the trajectory propagator. 

These flags are used to select the trajectory propagation 
scheme. A variety of schemes are available. See 
reference 2 for a description of the schemes. 


*1039-1040 


1044 


1079 


These flags define the coordinate system of the state 
output from the trajectory propagation portion of the 
program. 

This flag is used to select the desired MAESTRO analysis 
mode. Descriptions of the modes available are presented 
in the following sections. 

Flag used to read initial state from GTDS 24-hour hold file. 


1087 Flag used to read initial attitude from Attitude Determination 

Program. 

1047 Flag used to key doppler analysis. This analysis presents 

histories of thrust, weight acceleration due to thrust^ and 
doppler shift during a motor bum. 


* Flag must be input. 
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3 . 7 Optional Startup for Methods 5 and 8 

The numerical averaging methods can be started with osculating orbital elements 
by setting the value of KAVST (location 1098) equal to 1. Normally, it is assumed 
that the input orbital elements for these methods represent an average state. 

If KAVST is set to 1, however, the program performs a one -revolution integration 
of the osculating elements to obtain a set of mean elements for starting the 
averaged equations integration. The time corresponding to these mean initial 
elements is the average of the time for the first revolution or r/2, where 7 is 
the orbital period. 

The initial time is advanced from zero to 7/2, and the averaging integration is 
started with the mean elements at this new epoch. Print times and compute steps 
are advanced by the same amount and the output is given at times one-half revolution 
later than they would have if the averaged equations integration had been started 
with input mean elements. 

3.8 Shadow Calculations 

The shadow times are determined in two different manners. In one method, 
the time of shadow crossing is determined from an iteration procedure in the 
numerical integrator. The distance from the shadow cone is determined at each 
compute step. Tests are made to determine of the shadow cone was crossed or 
a minimum to the cone passed on the last step. If one of these conditions is 
satisfied, an iteration is made to find the time of crossing. The iteration consists 
of a Newton -Raphs on hunt to find the time which results in the spacecraft’s 
position at the shadow cone. During this hunt, the spacecraft's position is 
determined by interpolation. This shadow scheme can not be used with the 
averaging methods or multiconic. This method is also used to determine occultation 
times. The flag in location 1034 is used to select this option. The flag is set to the 
planet number where the observer is located. 

The other shadow method uses the osculating orbit to determine the intersection 
of the orbit with the shadow cone. The method uses the assumption that the 
osculating orbit will not change during one revolution. This method can be used 
with all of the trajectory propagation schemes; however, the orbit must neither 
impact nor escape the planet. 
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The flag in input location 1049 is used to select the desired shadow calculation 
technique. 

3.8 Run Stopping Criteria 

The run is stopped either by the time reaching the final time, location 4, or 
the spacecraft flying through closest approach to , or farthest distance from, 
the target planet. The flag in location 1032 indicates that a closest approach 
or farthest distance determination is desired. The 1079 flag determines whether 
it is to be closest approach or the farthest distance. The test for closest approach 
will not be attempted until the time since state epoch is greater than the time 
input in location 513. 

3. 10 Sample Case 

A sample run using the trajectory propagator is shown in Figure 3.2. The first 
part of the figure is a listing of the input data cards. The sample shown uses 
position and velocity vectors as the initial spacecraft state loaded into locations 
40 - 45, The epoch of the state is loaded into locations 50 - 55. The trajectory 
propagation scheme , used was the NICE /MEAN method and the automatic compute 
interval option was specified by the inputs in locations 1-3. 

The second part of this figure is the sample case output. The first part of the 
output presents a listing of the data cards read. This is followed by a complete 
listing of the input array including all of the program preset inputs. The trajectory 
output begins after the listing of the input array. This output consists of position 
and velocity vectors and magnitude in the designated coordinate systems. The 
orbital elements printed are semi-major axis, eccentricity, true anomaly, argu- 
ment of perigee, inclination and longitude of the ascending node. The output 
frequency flag in location 1030 is preset to zero, thus only first and last points of 
the trajectory were printed. The run stopped at closest approach to the target 
planet (Moon) because the KCRASH flag is preset to 2. 
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Table 3 . 1 

Trajectory Propagation Mode Inputs 
MODE = 0 

INITIALIZED INITIALIZED 


LOCATION 

VALUE 

DESCRIPTION LOCATION 

VALUE 

! DESCRIPTION 

I 

0 

error control limit 

1030 

0 

output frequency flag 

2 

30 

initial compute interval 

1031 

11 

target planet number ■ 

3 

30 

minimum compute interval 

1032 

2 

closest approach test fl* 

4 

540000 

Run stop time 



■ 

*30-35 

- 

initial s/c orbital elements 

1035 

0 

gravitational model nupaPr 

or 


or position and velocity 



40-45 


vectors 

*1036-1038 

4 

trajectory propagation 

38 

333.40 

initial s/c mass, kg 



method 

39 

0 

ephemeris time correction 

1039-1040 

3,2 

output coordinate systems 

47,48 

- 

s/c attitude, right 

. 1047 

0 

thrusting mode 



ascension and declination 

1049 

0 

shadow flag 

*50-55 

- 

epoch of s/c state 

1058 

0 

auxiliary output flag 

170-189 

n 

compute interval table 

1059 

0 

type of auxiliary output 

197 

4.7(10) 

solar pressure at 1 au 

1068 

6 

number of ordinates 

198 

.2 

s/c reflectivity 

1069 

3 

used in averaging 
number of intervals 

200-380 

- 

s/c motor characteristics 



used in averaging 

430-432 

- 

engine impulsive velocities 

1079 

0 

farthest approach flag 

433 

13560 

solar pressure area 

1084 

0 

solar pressure flag 

450-469 

- 

print table 

1087 

0 

input attitude element 
set number 

513 

0 

time to begin closest 
approach testing 

1088 

12 

unit number of attitude 
data 

*1001-1012 

S,E,M 

Planets in system 

1089 1234567 

satellite ID number 

*1014 

3 

numerical integration 
scheme 

1090 

0 

director’s display write 

*1015 

3 

launch planet number 

1094-1095 

0,0 

maximum number of jl 
zonals and tesserals us* 

*1017 

5 

ephemeris flag 



in gravity model 

1018 

1 

Earth oblateness flag 




*1019 

1 

input coordinate system 




1029 

0 

gravitational field flag 





* input required 
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.0 

0 .9 

0.9 

0.0 


0.0 

0 

.0 

0.0 

n ^ 

^ • . 


C(1 

J)= 0.10 4 43-0 4 

0 

.0 

0 .0 

0.0 

* 0.0 


0.0 

0 

.0 

0.0 

o 


C ( 15 , 

J)= -0 .2977D-Q 4 

0 

.0 

0.0 

‘ » 

0.9 


0.0 

0 

.0 

0.9 

0.0 


S{ 3, 

J)= 0.0 

-0 

.72130-05 

9.45370-05 



S( 7, 

J) = 0.0 

0 

. 14210-05 

0,5743^-96 

-9 .29190-05 


S( 4; 

J)* 0.0 

0 

.32990-05 

-0.23890-05 

-9.52230-96 

0 . 42480-95 

S( 5, 

J) = 0.0 

-0 

.69250-05 

-0 .317 80-26 

0.11590-95 

0.6296D-07 


-0 . 30 1 2D-0 7 






S{ 5, 

•J)= 0.0 

0 

.53420-05 

-n . H57D-Q5 

-9 . .34 3 40-06 

-0 .98790-07 


-0 .16 54D-07 

0 

.23730-03 




S( 7, 

J)= 0.0 

-0 

. 15790-05 

2 . 1 15 2D- 06 

0 . 24400-06 

0 . 1795D-07 


0 . 55 82D-0 3 

0 

.12270-08 

-0 , 91 3 ID- 10 



S( 8, 

J)= 0.0 

0 

.62070—05 

-0 . 1767D-05 

-0.14370-06 

-9 , 1757D-07 


0.239 1D-0 8 

-0 

.73900-09 

-0 . 5 9.13D-10 

0 , 979 2 D - 11 



29 . 



O L' 10 CO 


.■XSRS FIELD-POST MARINE 

C( 2, J)= -0 . 19 6 0 D-0 2 
( ?, J)= -0.3050D-OA 
( &, J) = -0.? 1*0 D-0 d 
( 2 , J) = 0.0 

( 0; J) = 0 ,0 


: 19 71 , I'! MOD 

= A , MMOD 

0.0 

■0 . 5400n-0<4 

0 .2 770 D-0 5 

0.0 

0.0 

O.o 

0.0 

0 .7900D-OA 

- m - v» - • • » i/ 1 ^ 

A A 

V . V 



Table 3.3 


Preset Physical Constants 


PLANET 

GRAVITATIONAL 

COEFFICIENT 

(KM 3 /SEC 2 ) 

MEAN 

RADIUS 

(KM) 

ROTATION 

RATE 

(RAD/SEC) 

SPHERE OF 
INFLUENCE 
(KM) 

Mercury 

2.168553xl0 4 

2.329816xl0 3 

0. 

0. 

Venus 

3.24853x10® 

6.098636xl0 3 

0. 

6. I6xl0 5 

Earth 

3.986008xl0 5 

6.37814xl0 3 

7. 29212361xl0“ 5 

9.25x10 s 

Mars 

4. 2915515x1 0 4 

3. 40953x1 0 3 

7,0882 xlO' 5 

5. 65x10® 

Jupiter 

1.267106xl0 8 

6. 98389xl0 4 

1.7 58549x10 " 4 

4.8xl0 7 

Saturn 

3.791869x10^ 

5. 7505x1 0 4 

0. 

5.4xl0 7 

Uranus 

5.80329xl0 6 

2. 5484xl0 4 

0. 

5. lxlO 7 

Neptune 

7.02607x10® 

2. 4983xl0 4 

0. 

8. 6xl0 7 

Pluto 

3. 31788xl0 5 

6. 3450x1 0 3 

0. 

3. 33xl0 7 

Sun 

1.3271545X10 11 

7.06xl0 5 

o. 

l.xlO 10 

Moon 

4.902778xl0 3 

1.73809xl0 3 

2. 6616995x10'® 

1. 1x10® 

Oddball 

0. 

0. 

0. 

0. 
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Table 3.3 (CONT’D) 


STATION NAME 

JOHANNESBURG 

TANANARIVE 

CARNARVON 

ORRORAL 

FAIRBANKS 

ROSMAN 

GREENBELT 

SANTIAGO 


* geocentric 


Preset Observation Stations 
LONGITUDE* LATITUDE* 


27.707280 

-25.735782 

47.303168 

1 c onooo/i 

-XWi U OUMUt 

113.716372 

-24.757607 

148.956961 

-35.448009 

-147.511732 

64.823643 

-82.876093 

35.013815 

-76.842965 

38.810366 

-70.666519 

-22.975938 
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input data cards 


TRAJECTORY PROPAGATION SAMPLE CASE 


50 6. 

53 15. 

40 10 51,8313 
45-10. 127676 
1 .0001 
10 56 4. 


51 10. 

5/5 1/5. 

41 50 32.0 a 8 3 
44 1.0265591 


52 1973. 

53 43.240 

42 2604.2559 
45-0,, R 0 7521 o 


2 30. 


2 20 . 


maestro output 


DATA CARDS FOR CASE l 


20 

0.600000000D 

01 

21 

o.ionoonoooD 

02 

22 

0.197300000D 

0 4 

23 

0 .150000000D 

02 

24 

O.lOOOOOOOOn 

01 

25 

0. 1250000003 

01 

50 

0.600000000D 

0 1 

51 

0 . 100000 0 00 D 

0 2 


0 . 1? 7 1000 OOP 

04 

53 

0 . 150000000 D 

02 

54 

0 . 1 ^ 0 0 0 0 0 0 0 D 

02 

55 

0 .4326°0009D 

02 

40 

0.10519.11100 

0 4 

41 

0 . 59 120/5 R 70D 

04 

47 

0 .7 6047559 OD 

04 

A 3 

-0.10 12 76 760 n 

02 

44 

0 . 3-026 5 5? ion 

01 

4 5 

-0 ,2 20 7 5 ?1 ROD 

01 

1 

0 . loooooooon- 

■0 3 

2 

0. 100000000D 

02 

3 

0 .20^0000000 

02 

10 16 

0 . 400 000000 D 

01 

0 

0.0 


0 

0.0 


LUNAR 

FIELD , NMOD= 0 MMOD=*** 






INITIAL JULIAN DATE 2441944.1352 
INPUT COMMON 



0 . 100 00000 OD- 

-0 3 

2 

0 .IOOOOOOOOD 

02 

1 

0. IOOOOOOOOD 

0.7 

A 

0.540009000D 

06 

5 

0 . IOOOOOOOOD 

01 

5 

0 . IOOOOOOOOD 

01 

K 

0 . lo 00 00 90 on- in 

10 

0 . IOOOOOOOOD 

21 

20 

0 .^OOooOOOOn 

01 

21 

0 . loooooooon 

02 

22 

0 . 19730O000D 

04 

? 7 

9.15000O000D 

02 

2 A 

0 .190000090D 

01 

25 

0 . 1 ? 5 r« 0 0 9 0 0 D 

01 

97 

0.4TR4417R7D 

05 

3R 

0.333 39 0099D 

03 

40 

0.1061R3130D 

04 

41 

0,6977048. ROD 

04 

42 

0.260425 5 90D 

04 

43 

-0 . 17*12 7 6 760 D 

0? 

AA 

0.70 76 559 ?OD 

01 

45 

-0 .2 80 35 21 ROD 

01 

46 

0 * 1 1744135 Of) 

05 

An 

0.167 161.7R8D 

0 3 

4R 

-0.14 854 429 3D 

02 

50 

0 . *oooooo^on 

Oi 

53 

0. IOOOOOOOOD 

02 

52 

0. 197300000D 

04 

5 3 

0 . .150000000D 

02 

54 

0 .140000000D 

02 

55 

0. 4326800900 

02 

10 0 

0.216P55300D 

OI 

101 

0 . 72 4 R 1.5 40 OD 

06 

10 2 

0.39 R60 3200 D 

06 

10 3 

0 . 429 155 150 D 

05 

10 4 

0 . 126710600D 

09 

105 

0. 379 1 R6000 D 

OR 

109 

0. 13271 5450D 

12 

110 

0 . 4.90777 7900 

0* 

112 

0.2 329 8 1600 D 

04 

113 

0 .S09863P00P 

04 

114 

0.677816 50 ID 

0 4 

115 

0.3409 5 3000D 

04 

116 

0 . 71^77D000D 

05 

117 

0 . 575O5OOO0D 

05 

HR 

0 .254R40000D 

05 

119 

0 . 249 93000 OD 

05 

ilf) 

0 .6R4-500000D 

04 

121 

0.706000000D 

06 

122 

0 . 17RO00000D 

04 

176 

0. 7 2923 L2 36 l.D- 

•0 4 

12 7 

0.70 R820000D- 

-0 4 

12 R 

0. 175 85 490 0 D- 

-0 3 

13 a 

9. 2663699 50 D- 

-05 

170 

0.36000 0000 D 

04 

171 

0.180000000D 

05 

170 

0. IOOOOOOOOD 

20 

180 

0 . 300000000D 

0 3 

181 

O.IROOOOOOOD 

04 

137 

O.IROOOOOOOD 

05 

19 7 

0.470000000D- 

-0 4 

1 QR 

0 .210000000D 

00 

70 1 

0 ."00OO00O0D 

OR 

20? 

0.200000 9 00 D 

05 

221 

0, 190000000- 

02 

2?2 

0.2000000000102 

260 

0.5 3 37 84000 D 

0? 

261 

0.266P92000D 

02 

7^7 

0 . 766 89 9000 D 

0? 

280 

-0 .9 78500 00 OD 

0 4 

2 8.1 

-0.97R500000D 

04 

30 1 

0 .900000000D 

0 3 

322 

0 .20 00000 00D 

0 5 

331 

0. 190000000D 

02 


0.7000000003 

02 

350 

0 .240400COOD- 

■01 

351 

0 . 1 20 7oo 009 n- 

-01 


0.i' , 9'>00OO0D- 

■01 

Sou 

u . 3532utu’UUD 

Ul 

361 

0 . 35 32000 OOP 

01 

390 

0 . IOOOOOOOOD 

21 

381 

0. IOOOOOOOOD 

21 

382 

0. IOOOOOOOOD 

21 

393 

0 .IOOOOOOOOD 

21 

384 

0 . innnnnooOD 

21 

385 

0. IOOOOOOOOD 

21 

£0 7 

O.IROOOOOOOD 

00 



40 P 

n T ^ 

09 

4 11 

0 . °3 55960 7°o 

00 

414 

0 . 3*79 36 19 9 Or 

01 

4l*» 

0.5049371 9 Qi 

01 

477 

0. 39-6009 9 00 D 

06 

4 34 

0 . 7700090009’ 

04 

4 3 7 

O.-’OOOOOOOOn. 

-01. 

440 

O. 7592000000 

06 

44 7 

0 .705 79 7900 P 

9 ? 

446 

9. 116500000D 

02 

45 n 

9. 1990000090 

21 

4 7 1 

n 1 £ "7 01 1 rt rt rt r\ 

' ♦ - ♦ ■ - ' r .f ^ KJ ■.» 1/ U 

rtrt 
U V.' 

47/ 

9.5909900000 

O 1 

479 

0 . 300999000D- 

-0 7 

4-33 

-0 .4 3.7101 75 7p 

00 

4Q^ 

0. 61 1 1.96 3550 

09 

490 

0 . 1000900990 

02 

4 o 7 

9 . 1.900999090-0 3 

4Q6 

0 .5909000000 

01 

59'’ 

9 .O^ 5 9 099900 

96 

50 5 

0.5 490000090 

03 

509 

n . 3300000990 

0 9 

190 " 

0 . 399907900D 

01. 

19 1 4 

0 . 7090 09 999 r> 

9 1 

19 1 p 

0. 1090909900 

01 

10 3 1 

n. 1109000000 

0.7 

19 3 C 

1 . 40001070 on 

01 

10 A 1 

9 . 1000 0 9 90 9 n 

0*1 

19 4 9 

9 . 199900999D 

01 

195 3 

9 . 590909900D 

n p 

1060 

0 . 19999999010 

01 

106 3 

9.709099900D 

01 

1066 

0 . 1090090000 

01 

1959 

9 . 399000990D 

01 

19 ’’9 

0.600009909D 

01 

10 9Q 

0. 133456 700D 

07 


AS9 

9 . 409009100.0 

09 

A 1_? 

0 . 1.9 34725 1 9n 

01 

A 15 

9.43357?4e in 

OT 


0.7070000030 

14 


1 . *799.799001 

00 

425 

9 . 709 900 099 n 

00 

ass 

9 . 799 30 01 OOP- 

-03 

4 A 1 

0.926090909P 

07 

44 4 

0 .293999009^ 

0 4 

447 

0 . 6090000090 

04 

At 0 

0.i9i099999n 

95 


0 .2040 DO 090 P 

0? 

475 

0.71O00O09OP 

12. 

439 

-9 . 449 l 74 1 31p 

99 

433 

-0 .61.06 944 7017 

09 

4 36 

0 . 6 17 06 OC 7 1 ri 

00 

401 

0 . 1091000900 

9 3 

40 4 

9 . 199900191-)- 

-0 3 

407 

0.2 r '9000971n 

00 

59 7 

9 .565900001P 

06 

506 

o.5i 00000900 

0 0 

509 

9 . 1 99900990P 

1.1. 

10 10 

9.199999059^ 

1 1 

19 13 

0 . 7599999011 

0 1 

1.0 19 

9 . 1990099900! 

01 

10 32 

0 .2010000901 

01 

10 79 

0 . uniOOOOOn 

01 

10 43 

9 . i 599903900 

07 

10 5 1 

9 . 109009 9 OOP 

0 2 

1054 

9.1700000000 

0 2 

1961 

3 . 11990090019 

07 

196 4 

0.1990999901 

01 

19S* 7 

9. 19 9900 910 p 

0 3 

19 70 

0 .9509003000 

02 

10 3? 

0. 11.909OOOOP 

07 


A 10 

1 . * 3 76 OHS 71 

90 

41? 

9 . 76901904] p 

01 

A 1 £ 

9 . /0/7077661 

n -1 

A "''* 

9 . 1 165110011 

97 

A ?? 

9.1 7569179Q1 

03 

ASS 

0 . 7090990001- 

-91 


0 . 7600000901 

95 

4 A 2 

0,73730 9 noop 

07 

A A S 

0 ,1.31 4361 76p 

01 

£48 

79.600000000P 

04 

4 70 

J . 61990006 JD 

01 

47? 

9 . 173.37710OP 

07 

47R 

•9 . 770000QQ0 D 

0 4 

AS 1 

-9 . 771173499^ 

09 

A OA 

9 . 1 !. 31.9 96090 

01 

A SI 

-1 .5765706971 

00 

A Q 7 

0 . 1 009999000 

07 

495 

9.79990OO9OP- 

-91 

571 

0 „ * 1 60 90091 1 

06 


9.4309199909 

00 

507 

9.9601001001 

03 

5 1 0 

9 . 11.009 09 OOP 

06 

10 11 

9 . 799990909^ 

01 

10 1 n 

n .301099900-) 

01 

19^3 

9. 10099900 on 

01 

10 ^ 5 

0 . 1909099001 

01 

10^0 

9 . 7990909991 

01. 

10 A 3 

0. 719090 0 901 

02 

10 5° 

9 .2009000000 

01 

10 5 5 

0 . 1 O09O9H9P 

01 

10*0 

9.70O9OOO09D 

0.1 

10 5 5 

9 . 1.99O09090P 0? 

10^0 

9.6909990900 

01 

10 71 

0 . 6900090900 

01 

10 Q 8 

0 , 1 709O99OOD 

07 
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time= 

earth 


0.0 


MEAN EQUATOR AMD EQUINOX OF 1.950 
0 • 105 1R313D 04 0.59130403D 04 0.?604?559n 04 0.654336000 04 

-0.10127676D 02 0 . 30 25 550 ID 01 -0,780152190 01 0.1093S704n Q7 

ORBITAL ELEMENTS W.R.T. EARTH. 

0.2.12754400 06 0.969150530 00 -0 . 60 7953600-07 0.127966500 07 

0.2 9190 3 4 on 0 2 -0 . 4529 a 14 80 0 2 

MOON TRUE EQUATOR AND PRIME MERIDIAN 

0.39 4 926910 06 0.493004-790 05 0.457519780 05 0.400392900 05 

-0.92660062D 01 0.528730230 01 -0. ’40 6.38106 0 01 0 . 1 . 1 41.6 177u 02 

ORBITAL ELEMENTS W.R.T. MOON 

-0.375255440 02 0.65905917.0 04 -0.517440500 02 0 . 23 6 ?7 D 03 

0 . 2641.4 79 Ro 02 -0 . 1596.36640 0 3 

0 4 

CLOSEST APPROACH TO TARGET PLANET 
TIME- 110.0075 

EARTH MEAN EQUATOR AND EQUINOX OF 19 *0 

-0 .95471 .3640 05 -0,. 35 86 3. 126 0 06 -0 . 16 90 77 76r> 06 0 . 4 f) 740 50 on O* 

-0.379328120 00 0.911527410 00 -0.102784970 01 0.14216765D 01 

pRBITAL ELEMENTS W.R.T. EARTH 

-0.62071096^ 07 0 . 106395 LID 01 -0.229564730 02 0 . 7?9620 4 »r>. 02 

0.113173510 07 0.862 70 6 05 D 0? 

MOON TRUE EQUATOR AND PRIME MERIDIAN , 

0.732540560 00 0.1405416?n 01 -0 . .14].6 «5?4n 01 0.207717840 01 

ORBITAL ELEMENTS W.R.T. MOON 

-0 . 71 359554D 04 0 . 13555.1 iAr> 0.1 0.173077700-01 . 0.1411 4. 5 mu 07 

0.11594.3370 0 3 -0.7479666 90 02 

36 7 4 

NO MORE DATA, RUM TERMINATED 
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Section 4 

RETRO- FIRE ANALYSIS MODE 

This mode is used to determine the retro motor's firing time and attitude. When 
the MODE flag is set to one, the entire analysis is performed using numerical 
integration. MODE set to 2 option uses some conic approximations to propagate 
the state. See reference 1 for a detailed description of the approximations. 

The retro motor's firing time and attitude are determined by making repeated trial 
retro firings at various times and attitudes. The firing times and attitudes are 
scanned in a systematic order so that displays can be presented. The user can. 
then select the desired firing time and attitude by studying the displays. The 
analysis mode does not automatically select a firing time or attitude. 

The user has three available options to scan the firing time and attitude space. 

The desired option is selected through the KAPOPT flag as follows: 

KAPOPT = 1. Firings are made between an input range of true anomaly 
about perigee on the approach hyperbola. The number of 
firings is determined by the KFIRE input in location 1043, 

The attitude range is an input cone angle about the initial 
right ascension and declination in locations 403 and 404. If 
these values are zero, the cone is centered about the velocity 
vector at the final time. 

2. Firings are made from asymptote to asymptote on the approach 
hyperbola. The number of firings is selected by the KFIRE 
input in location 1043, The attitude range is input through 
initial right ascension and declination, increment in right 
ascension and declination, and number of right ascensions 

and declinations to be tried. These values are input via 
locations 403 to 406 and 1041 to 1042. If locations 403 and 404 are 
both zero, the initial right ascension and declination is defined by 
the velocity vector at closest approach to the target planet. 

3. Firings are made between input times from liftoff. The firings 
are made between the initial time in location 400 and the last 
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time in location 401. When MODE = 2, the range is divided 
into KFIRE steps of constant true anomaly in the range. When 
MODE = l,the range is divided into KFIRE number of equal 
time steps. The attitude range is input the same as described 
in 2. 

If MODE = 1 option is used, KAPOPT = 3 must also be used. 

When the retro motor is to be numerically integrated (MODE = l),its thrust and 
weight characteristics should be loaded into engine number two. 

The RAE-B retro-motor characteristics are pre-set in engine 2. Table 4.1 
presents these characteristics. 

This mode has the ability to save the post-retro state on an auxiliary unit for use 
in other analyses. The flag in location 1082 controls this option. It is set to 
the desired unit number. 

If the flag in location 1092 is nonzero, the post-retro state for every firing will be written 
on an auxiliary unit using subroutine PUTELS for retrieval by the GTDS program. 

The inputs preset in the trajectory Propagation Mode also must apply in this mode. 

Table 4.2 presents the necessary inputs for an approach analysis plus the preset 
value of the inputs pertinent to this mode. Notice that it is possible to obtain the 
initial state from the previous Midcourse Correction Analysis. This is accomplished 
by setting KREAD (location 1057) to the midcourse number of the desired execution 
time. A sample approach analysis is shown in Figure 4.1. The first part of this 
figure contains the input data. For tliis rim, it consists of the initial date, MODE 
set to 2 for an approach analysis ,and the KREAD flag set to 2 to pick the initial 
conditions at the second midcourse execution time. The second part of this figure 
is the approach analysis output. A listing of the input is presented first. This 
is followed by the input array with the preset values. The approach analysis 
output shown is the minimum eccentricity grid and minimum trim fuel grid. The 
value in the first grid is the minimum lunar orbit eccentricity for all the firings 
made at the indicated attitude. The value in the second grid is the minimum trim 
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fuel required to obtain the desired lunar orbit for all firings made at the indicated 
attitude. The desired lunar orbit is defined as a circular orbit at the radius and 
inc lin ation with respect to lunar equator described in locations 444 and 446. Thus 
a quick scan of this grid should indicate a range of desired attitudes. 

A second sample of the approach analysis is shown in Figure 4.2. This run is 
similar to the last except that the KAPOUT flag (1046) is set to one and fewer 
attitudes are studied. The KAPOUT flag set to one will print information about 
each firing on the approach hyperbolas. The fewer attitudes were used to reduce 
the output. If a 15 x 10 grid were used, as in the first run, 150 pages would be 
printed. The output from this run is similar to the other runs in that both the 
input cards and input array is presented. The next part is a firing time analysis 
at the attitude shown. This output presents information about the lunar orbit for 
each firing time. This output is divided into two parts for easier readability. 

The first part contains 


FIRING 

TIME 

Time the retro is fired, time since liftoff 

NEXT 

SHADOW 

Time until shadow. If this quantity is neg- 
ative, it is the time until there is no shadow. 

SHAD FREE 
TIME 

Time between shadows 

DELTA V 
TRIM 

Trim maneuver to correct the inclination to 
the desired value and circularize at the desired 
radius. 

TRIM 

FUEL 

Fuel used in the trim maneuver. 

NORMAL VECTOR 
RT, ASC DEC 

The right ascension and declination of a 
vector normal to the orbit plane in mean 
Earth equator and equinox of 1950. 


The second part contains the elements of the post-retro lunar state with respect to 
the true lunar equator and node at the time of closest approach. 


FIRING 

TIME 


Time the retro is fired, time since liftoff 


SEM 

ECC 


Semi-major axis of the lunar orbit. 
Eccentricity of the lunar orbit. 



TRUE 

ANOMALY 

AOP 

INC LIN 

LON 


True anomaly of the spacecraft in the orbit. 

Argument of pericynthion 

Inclination of the lunar orbit to the lunar 
equator. 

Longitude of the ascending node. 
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Table 4. 1 


Retro Motor Thrust/Weight Characteristics 


Time since 

Thrust 

Mass flow 

ignition (sec) 

(Newtons) 

rate (KG/SEC) 

0 

-9785.0 

3.532 

19 

-9785.0 

3.532 

20 

0 

0 
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Table 4,2 

* 

Retro Fire Analysis Mode 
MODE = 1 or 2 


INITIALIZED 


LOCATION 

VALUE 

DESCRIPTION 

20-25 

- 

liftoff epoch 

**30-35, 40- 
or 1057 

-45 - 

initial state 

220-386 

- 

retro motor in engine 
2 for MODE = 1 

400 

- 

initial firing time 

401 

- 

final firing time 

402 

- 

increment between firings 

403 

- 

initial right ascension 

404 

- 

initial declination 

**405 

- 

increment in right 
ascension 

**406 

- 

increment in declination 

407 

.68 

impulsive velocity of 
retro 

410-419 

- 

observation site longitude 

**441 

226 

midcourse motor ISP 

**442 

282.5 

retro motor ISP 

**443 

71.11 

mass of retro fuel 

**444 

2838 

desired final orbit radius 

k*446 

i 

t 

59.0 

desired final orbit 
inclination 

473 

13.77 

retro drop weight 

474 

5 

attitude cone angle 

|475 

20 

true anomaly firing range 

480-489 

- 

observation site latitude 


I 

i 

\ 

\ 

* inputs for trajectory propagation also set 
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LOCATION 

VALUE 

DESCRIPTION 
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0 

retro motor burn time 
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2 

retro motor number 

**1041 

10 

number of right ascensions 
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15 

number of declinations 

**1043 

20 
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- 
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1 
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1 
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11 
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- 
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1081 

0 
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0 
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ysis write set to unit no, 
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- 

Element set no. of input 
attitude data 

1088 

12 

attitude unit number 

1089 1234567 

satellite ID number 

1092 

- 

Element set no. of 
write for GTDS 
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FIGURE 4 . 1 


PJ’T P.0 

minimum 


MOTOR DETERMINATION SAMPLE CASE 
ECCENTRICITY AND MIMIVtt' mEL CHIP'S 
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ATTITUDE SCA> T 


mini: :um eccentricity 


AN CHO 75 VECTOR EPOCH 
IS. HR lit. MTm 47.36QSHC 
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R. 
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Figures *.2 


INPUT CARDS 
20 6 . 
22 15. 
50 6. 

5 2 15. 
10*1 1 . 
1057 1. 


MAESTRO ANALYSIS 


DATA CARDS FOR CASE 
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FIRING 

TIME 

(HR) 

109 .990 

109.9 73 

109 .959 
10 9 .95 9 
109.990 
109.99 2 
110 .004 
110.017 
110.029 
110 .0*1 
110.053 
110 .055 
110.0 77 
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110 . 10 1 

110 . 113 
110 .126 
110.130 
110.151 
1 10 . 16 3 
110 .176 


FIRING 

TIME 
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110.017 
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110.0®9 
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110 .139 
110.151 
110.16 3 


APPROXIMATE APPROACH ANALYSIS 

ANCHOR VECTOR EPOCH 
15. IIR 14. MIN *3 . 26 °SEC 


DEC. R.A. 

-30.071 135.065 


NEXT 
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Section 5 


MIDCOURSE ANALYSIS MODE 

The midcourse analysis mode is designed to present to the user displays from which 
he can select the ,T best” correction according to his needs and constraints. The 
core of the midcourse analysis mode is a group of subroutines used to determine 
the correction at a specified place on the transfer trajectory which satisfies input 
end conditions, see reference 1. Logic was built around those subroutines so 
that certain control variables can be varied in order to automatically scan the 
control variable space. 

The KTF flag is used to determine the type of scan, and hence, the type of displays 
to be printed. This flag is defined as follows: 

KTF =0 A one -dimensional scan of midcourse 

execution time. The initial time is loaded 
into 478, the increment in 434, and the number 
of execution times in 1051. 

>0 Two-dimensional scan of midcourse execution 

times and flight times. Execution scan is 
defined as above. The flight times are scanned 
in KTF one-hour steps beginning at the 
desired flight time in location 422, Fixed time 
of arrival guidance must also be specified. 

< 0 Fixed attitude scan. A two-dimensional 

scan of midcourse execution times and mid- 
course impulsive velocities at a fixed attitude 
input via locations 47 and 48. The impulsive 
velocity is varied about the velocity input in 
location 426, in -KTF steps of size equal to the 
value in location 479. The attitude can be sys- 
tematically varied about the input attitude if desired 
The number of attitude steps are input in 1041 
and 1042. The attitude range is input in 474. 

When this logic is used, the attitude is varied 
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around a cone whose half angle Is input 
in location 474. The cone is centered at 
the input attitude. 

Besides the scan option listed above, the user also has the option of selecting 
one of five guidance laws. The user selects the desired law through the KGLAW 
input (location 1063) as follows: 


KGLAW = 1 

Minimum fuel. This is a critical plane 
maneuver correcting only B*T and B*R (*), 

KGLAW = 2 

Fixed time of arrival. This corrects B*T, 

B • R (*) and time of flight to the desired values 

KGLAW = 3 

Fixed Target Energy, This corrects B • T, 
B* R (*) and energy on the mean passage 
hyperbola to the desired values. 

KGLAW = 4 

Variable Target Energy. This corrects 


B • T, B* R (*) and the energy in lunar orbit 
to a desired value. The lunar orbit is 
determined by firing a retro motor at the 
perigee point of the passage hyperbola along 
the velocity vector. 

KGLAW = 5 Minimum Total Fuel. This corrects B * T 

B* R (*) while minimizing the total midcourse 
and trim fuel expenditure by optimizing retro 
orientation and firing time. 

The desired end conditions can be adjusted if a retro overburn exists. An overburn 
occurs when the available retro velocity is greater than the velocity to circularize 
at the desired radius. When this condition exists, it is necessary to use up the 
extra velocity in some way to avoid a post-retro maneuver. Two methods are 
available to the user. They are called the variable inclination procedure and the 
variable periapsis procedure. These procedures are described in detail in reference 
1. Simply, the variable inclination method adjusts the aiming point so that an 
out-of-plane retro burn at periapsis of the approach hyperbola will result in a 
circular orbit at the desired radius. In the variable periapsis procedure, the 
periapsis radius is adjusted so that an in-plane retro firing at the desired radius 
will result in a circular orbit at the desired radius. The IVTI flag is used to 
select either of these procedures as follows: 

(*) or, alternatively, radius of closest approach and inclination 
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IVTI =* 0 


in-plane retro maneuver antiparallel to 
velocity vector at periapsis. 

= ± 1 Variable inclination procedure approaching 

above or below desired inclination. 

= t 2 Variable periapsis procedure circularizing 

in-plane before or after periapsis. 

It should be noted that none of these options are executed when the NORMIN flag 
(location 1080) is non-zero. In this case, the retro-maneuver is determined from 
a steepest ascent optimization procedure. This procedure is explained in reference 

1 . 

The IBURN flag in location 1070 is used to select the method in which the midcourse 
motor is applied. The number of this flag corresponds to the trajectory propagation 
method defined in location 1036, Any of the methods may be used except the 
averaging techniques (input value of 5 or 8). It is suggested that Cowell be used 
when an accurate simulation is required. When less accurate simulations are 
required, multi-conic or impulsive velocity can be used. Impulsive velocity is 
specified by a zero setting of the flag. The midcourse motor's thrust and weight 
expended histories are input in engine number 1. The values for the RAE-B 
midcourse motor are pre -initialized in MAESTRO, These values are shown in 
Table 5.1. The impulsive velocity of the midcourse motor is determined from 
the standard rocket equation. The mass and specific impulse input via locations 
38 and 441 are used in this equation. 

The size of the midcourse velocity step in the midcourse determination is limited 
in order to help ensure convergence. The maximum step is the product of the 
limiting factor in location 1067 and the partial derivative step in location 479. 

Information about the midcourse correction at each execution time is stored on 
an auxiliary I/O unit. This unit number is input via location 1061, This input 
must correspond to a unit defined with the proper JCL. This information is 
retrieved in other analysis modes. 

The inputs preset in the Trajectory Propagation Mode are also preset for this 
mode, A guide for the use of this mode is shown in Table 5,2, This table 
presents the inputs pertinent to the Midcourse Analysis Mode, 


50 



A sample of the Midcourse Analysis is shown in Figure 5,1, The first part 
of this figure is the input data. For this run, the initial state is loaded as 
position and velocity vectors in locations 40-45. The initial date is in locations 
50-55. The KGLAW flag in 1063 is set for fired time of arrival guidance and 
the MODE flag in 1044 is set to 3 for Midcourse Analysis. The KTF flag 
(location 1077) is set to 2 to obtain a scan at two flight times. 

The first part of the program output, second part of Figui'e 5. 2, consists of a 

listing of the input cards followed by a listing of the input array including the 

preset values. The midcourse analysis output follows the listings of the input 

array, A line of information is printed for each execution time and flight time. 

Since it was necessary to abbreviate the titles, a description follows: 

MC Midcourse number 

NO 

TIME Execution time. Time since liftoff. 

MAG Magnitude of the midcourse correction. 

FUEL Fuel required for the midcourse maneuver. 

BURN Midcourse motor burn time. 

TIME 

R.A.. . Right ascension and declination of the correction 

DEC, velocity in the Earth mean equator and equinox of 1950. 

SPIN-SUN Spin axis-Sun angle. Angle between the Sun and the spin 
axis when the maneuver is applied. 

MIN SUN Minimum angle between the spin axis and the Sun during 
the reorientation to the midcourse attitude, 

ATT FUEL Fuel used in the attitude orientation maneuver. 

TRACKING The rounded off elevation angle for each tracking station. 
VISIBILITY If the elevation angle is between 0.0 and 10.0, a 1 is printed, 
etc. The station is identified by the first letter of its name. 

The Midcourse Analysis output also includes predicted errors in mission objectives 
due to midcourse maneuver execution errors. The answers are computed by 
linear propagation of ensemble statistics. The user supplies 
Location Name 

435 SIGAT lo pointing error in degrees (actually 

assumed to be 97% probability on Rayleigh 
distribution) . 
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436 StGDV lc proportional error in velocity (fraction 

of midcourse velocity in error). 

A 1c resolution error of . 3 meters/ second is built into the program in addition 
to the input proportional error. The user must also supply the probability 
(location 1070, in percent) to which he would like the output errors scaled, 
assuming a Gaussian distribution. The output format lists the following: 


MC 

Midcourse execution number. 

TIME 

Execution time. 

TFUEL 

Total fuel variation relative to errorless execution 
(primarily caused by trim requirements and assuming 
no later midcourse) 

PERILUNE 

closest approach distance error 

INCLIN 

Inclination error 

TFLITE 

Flight time error relative to errorless execution. 

VJNFIN 

Hyperbolic excess speed error relative to errorless 
execution. 

VCIRC 

Circular excess speed error relative to errorless 
execution. 

MC2 

Expected correction velocity requirement at a 
second midcourse. 


The second midcourse correction is presumed to use the same guidance law as 
the first. The time of the second maneuver is input via location 440 in seconds 
past liftoff epoch. 

Midcourse analysis output also includes lunar orbit insertion information computed 
from the predicted arrival conditions based on a successful midcourse maneuver 
execution. The listed output items are: 

MC NO Midcourse execution number 

RETRO TIME Best retro firing time based on minimum trim fuel 
RETRO RA Right ascension of the spin axis for the retro maneuver 
ATT DEC Declination of the spin axis for the retro maneuver . 
SPIN-SUN Spin axis-sun angle for the retro maneuver. 

TRIM VEL Trim velocity to achieve desired lunar orbit. 
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TRIM FUEL 

TOTAL FUEL 

TRACKING 

VISIBILITY 


Trim fuel 

Sum of the midcourse and trim fuel expenditures. 

Visibility for the retro maneuver — same format as 
for the midcourse maneuver. 


This sample case contains two lines of print at each execution time. The first 
line corresponds to the flight time loaded into location 422 while the second line 
has a flight time of one hour longer. These times can be verified from the retro 
firing times in the third block output. More than one flight time was used because 
the KTF flag in location 1077 was set to two. 


A second sample case using the fixed attitude option is shown in Figure 5.2. The 
data for this case is very similar to the last except for the special inputs pertaining 
to the fixed attitude mode. The KTF flag is set to minus three to specify 3 
velocity increments. The increment in velocity is input in location 479 and its 
mid-value input in location 426, The center of the attitude scan cone is the 
spacecraft attitude input via locations 47 and 48. The half cone angle is input 
via 474. The attitude cone is to be divided in three segments in both right ascension 
and declination as specified by locations 1041 and 1042. 

The MAESTRO output is in the second part of the figure. The input cards followed 

by the entire input array is presented first. The fixed attitude scan “follows. The 

fixed attitude scan output is divided into blocks pertaining to the midcourse motor 

firing time. In each block the following information is printed: 

DVM Midcourse impulsive velocity 

RTA right ascension of midcourse velocity 

DEC declination of midcourse velocity 

RCA radius of closest approach at target planet 

INC inclination of approach hyperbola w.r.t. target planet’s 

equator. 

TFLT time of flight, time since liftoff 

ROPA radius at the post-retro apsis opposite approach 

pericynthian, assuming a pericynthian maneuver 
opposing the arrival velocity. 

FCP Total fuel, assuming a trim maneuver which circularizes 

in-plane at approach pericynthian after pericynthian retro 
maneuver opposite to velocity vector. 
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TCF 


total fuel to trim to a circular orbit with the desired 
radius . No inclination change made. 


Notice that the output consists of a 4-dimensional scan of midcourse control 
variables - execution time, midcourse velocity, right ascension of midcourse 
velocity and declination of the velocity. 
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Table 5.1 


Midcourse Motor Thrust/Weight Characteristics 


Time since 
ignition (sec) 

Thrust 

(Newtons) 

Mass flow 
rate (KG/SEC) 

0 

45.5 

0.02404 

900 

26.6892 

‘ 0.01202 

20000 

26.6892 

0.01202 
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Table 5.2 

Midcourse Analysis Mode * 

MODE = 3 

INITIALIZED INITIALIZED 


LOCATION 

VALUE 

DESCRIPTION : 

LOCATION 

VALUE 

DESCRIPTION 

20-25 

- 

Launch epoch 

473' 

13.77 

retro drop weight 

**30-35 

— 

initial state 

474 

5 

attitude cone angle used 

or 





in fixed attitude node 

40-45 



475 

20 

true anomaly range used® 

200-386 

- 

mideourse motor in 



in retro determination ® 



engine 1 



when 1080^ 0 . 

407 

.68 

retro impulsive velocity 

**478 

7200. 

initial execution time 

410-419 

- 

observation site longitudes **479 

.0003 

velocity increment for 

**420 

2838. 

Desired lunar radius 



partials or fixed attitude® 
scan. ™ 

**421 

-59.0 

Desired lunar inclination 

480-489 

— 

observation site latitudes 

422 

396000. 

Desired flight time 

490 

10 . 

Tolerance on B»T ® 

423 

0.62 

Desired hyperbolic 
excess speed 

**491 

10 . 

Tolerance on B*R 

424 

0.0 

Desired circular 
excess speeds 

492 

10 . 

Tolerance on time of 
flight 

434 

7200. 

Increment in execution 
time 

493 

0.001 

Tolerance od hyperbolic 
excess of speed 

435 

.7 

Pointing error 

494 

0.001 

Tolerance on circular 
excess speed 

436 

.02 

Proportional error 

495 

.02 

Tolerance on minimum 

440 

259200. 

Second midcourse time 



total fuel 

**441 

226.0 

midcourse motor ISP 

496 

5. 

Tolerance on closest 

**442 

282.5 

retro motor ISP 



approach 

**443 

71.11 

Mass of retro fuel 

497 

.2 

Tolerance on inclination 

**444 

2838 

desired final orbit radius 

1030 

-1 

Output frequency flag 

**446 

59 

desired final orbit 
inclination 

**1036 

6 

Trajectory propagator 
scheme 

447 

6000 

Desired B • T 

1041 

10 

number of right ascension 
in fixed attitude scan 

448 

6000 

Desired B - R 

1042 

15 

number of declinations in 


fixed attitude scan 


* inputs for trajectory propagation also set 

** required inputs 
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Table 5.2 (Cont’d) 


INITIALIZED 


LOCATION 

VALUE 

DESCRIPTION 

1043 

20 

number of firings in retro determination, 
used when 1080^ 0 

1044 

- 

MODE flag set to 3 

1050 

0 

Auxiliary output flag 

1051 

10 

Nmnber of executions 

**1062 

2 

Miss vector option flag 

**1063 

1 

Guidance law (1. = min. fuel) 

1064 

1 

Number of trials to recompute secant matrix 

**1065 

10 

Number of trials allowed 

1066 

1 

Preliminary targeting done 

**1067 

100 

Limiting factor 

1070 

95 

Output probability 

1071 

6 

Integration type for midcourse burn 

**1075 

6 

Trajectory propagation method when generating partial derivatives 

1076 

- 

element set number of state from GTDS program. 

**1077 

0 

midcourse type flag 

**1078 

0 

overburn option key 

1080 

0 

retro optimization flag 

1089 

1234567 

satellite ID number 

1093 

- 

plot tape unit number 
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FIGURE 5.1 


MIDCOUP.SE ANALYSIS SAMPLE CASE 
EXECUTION AHD FLI^H" TIME SCAN 

INPUT DATA CARDS 


20 

6. 

2.1 

10. 

22 

1973. 

2 3 

15. 

24 

1 . 

25 

3.25 

50 

6 . 

51 

10. 

52 

19 73. 

5 3 

15. 

54 

14. 

55 

43.268 

40 

1051. 8313 

41 

5932.0483 

42 

2604.7559 

4 3- 

■10.141568 

44 

3.03071049 

45- 

•2.80 736 7 

104 4 

3. 

106 3 

2. 

1077 

9 

* • m 

1051 

5. 






MAESTRO OUTPUT 


DATA CARDS FOR CASE 1 


20 

0 .600000000D 

01 

21 

0. 100000000 D 

02 

22 

0 .397300000D 

04 

23 

0 . 150000000D 

02 

24 

0.100000000D 

01 

25 

0.3950000000 

03 

50 

0.600000000D 

01 

51 

0 . 100000000D 

02 

5? 

0 . 3 97 300 000 D 

04 

53 

0.150000000D 

02 

54 

O.74OOOO0OOD 

02 

55 

0 .4096800000 

0 ? 

40 

0 . 10 5183130 D 

04 

41 

0 .593204R30D 

0 4 

42 

9.9^049559OD 

04 

43 

-0. 1014156 80 D 

0 2 

44 

0 , 30 30 810 4 9 n 

0i 

45 

-0 .2P0706700D 

01 

10 44 

0 . 300000000D 

01 

106 3 

0 .7000000000 

01 

10 i 8 

0 .200000000D 

01 

10 51 

0.500000000D 

01 

0 

0.0 


0 

0.0 


LUNAR 

FIELD, NMOD= 0 

MMOD= 

0 






INITIAL JULIAN DATE 2441844.135? 
INPUT COMMON 


A 

0 .540000000D 

06 

5 

0.100000000D 

0.1 

6 

0 .1000008000 

01 

8 

0 . 1000000000-10 

10 

0 . 1000 00 GOOD 

21 

20 

0.6000000000 

07 

21 

0 .100000000D 

02 

22 

0 . 19 7 30 0000 D 

04 

23 

0 .1500000000 

02 

24 

0 .100000000D 

01 

25 

0.3 2 5000 GOOD 

01 

27 

0,4] 94 47 9 5 ->d 

05 

38 

0. 333390000D 

03 

40 

0.105133] 30 D 

04 

41 

0.59 9204 9500 

04 

42 

0 .2604255900 

04 

43 

-0.10 14 156 30 D 

02 

A A 

0 . 3O3091Q49D 

01 

45 

-0.280736700D 

01 

46 

0.4 7.944 13520 

05 

47 

0 .1633611980 

01 

48 

-0. 14 854 42 76 D 

02 

50 

0. 600000000 D 

01 

51 

0.700000000D 

02 

52 

0 . 19 7 30000 OD 

04 

53 

0. .150000000D 

02 

54 

0.1400000000 

02 

55 

0.432680000D 

02 

10 0 

0 .216 35 5300 D 

01 

101 

0 . ?9499.54Q0 D 

06 

10 2 

0. .3936032000 

06 

10 3 

0 .4291551500 

05 

1,2 A 

0. 11617.06000 

09 

105 

0.379186900D 

03 

10 9 

0 . ] 32715 450 D 

12 

110 

0 .4909717900 

04 

112 

0.2 3 298 16 00D 

04 

113 

0 .60936 36 00 D 

04 

114 

0 .63781.650 1 0 

04 

115 

0 . 34095 3000D 

04 

116 

0 .7142200000 

05 

1.17 

0 .5150500000 

05 

118 

0 .254840000D 

05 

119 

0 .2498300000 

05 

120 

8.6 3 45000 00D 

04 

121 

0 . 706000000 D 

06 

12 2 

0. 173800000D 

04 


0.1797113610- 

•04 

127 

0 . 708829000D- 

•0 4 

12 8 

0 . 175 3549000- 

•8 3 

134 

0 .765 1 <799590- 

■06 

170 

0 . 36 0000000 D 

06 

171 

o.ioooooooon 

20 

130 

3 .1S0000090D 

05 

X 8 1 

0.180000000D 

04 

19 7 

0.4700000000- 

■0 4 

19 3 

0.7 7 0000 GOOD 

00 

201 

0.900000000D 

0 3 

202 

0 .200000000D 

0 5 

9 9 \ 

0 .1900000000 

0? 

n " 

^ n 1 'i '3 * ' i 3 

S 

U Z, 

26 u 

u .5337o4uuuu 

u2 

2bi 

u . 2 bt> yy 2'jo ud 

U2 

262 

0.266892000D 

02 

2 80 

-0 .9785000000 

04 

2 81 

-0 .9 78500000 D 

04 

321 

0 ,9O9°0O00OD 

0 3 

322 

0.200000000D 

05 

331 

0 .190000000D 

0? 

332 

0 . 200000000D 

02 

350 

0 .2404000000- 

•01 

351 

0 .7.202000000- 

■01 

352 

0.120200000D- 

■01 

360 

0 .35 3200 GOOD 

01 

.361 

0.3532000000 

07 

380 

0 . 100000000D 

21 

381 

0 .100000000D 

21 

382 

0 .100000000D 

2 7 

383 

0 . 10000 0000 D 

21 

. 3 84 . 

. n .1000000000 

21 

385 

n inononnnnn. 

7J. 
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40 7 

0 .630000000D 

00 

40 8 

0 . 7 36090000D 

09 

409 

0 .490000000D 

09 

410 

0. 483583260D 

00 

AH 

0 . 32559602 3D 

00 


0 .19 34725 10 D 

01 

413 

0 .2 59 9 7 89 4 ID 

01 

414 

0 # 370^619900 

01 

415 

0 . 4 076 7 ? 46 1 D 

-91 

416 

0.494202256D 

01 

4X7 

0.50^9 S2138P 

01 

429 

0 .2 83Q09000D 

04 

421 

0 . 116500000 D 

0 3 

422 

0 . 396000000 D 

06 


0 .620000000D 

00 

4 33 

0. 1356 00 000 D 

05 

4 34 

0 .720000000D 

04 

435 

0 . 7OO90OQOOD 

00 

436 

0.200000000D- 

-01 

437 

0 . 7 0 0 0 0 0 9 0 0 D- 

-01 

4.3 8 

0 .3OOO8OO00D- 

-03 

439 

0.3600000 00 D 

05 

440 

0 .259 200000D 

06 

44 X 

0 .226000000D 

93 

442 

0.282500000D 

0 3 

44 3 

0.705 79 2000D 

02 

/[A A 

O.'iRiROOOOOD 

04 

445 

0 .131*361760 

01 

446 

0 . 116 5090 00 D 

03 

447 

0 .6090000000 

04 

448 

0 .600000000D 

04 

450 

0 . 1000 9 0009 D 

21 

460 

0 . !°OOO0OO0D 

0 5 

470 

0.600000000D 

01 

471 

0 .16*7 8000000 

00 

/. 70 

0.^040000000 

02 

473 

0 . 12 337700QD 

02 

474 

0 .500090000D 

0.1 

4 75 

O.IOOOOOOOOD 

02 

478 

0 .720000000D 

04 

4 79 

0 .300000000 D- 

-0 3 

4 80 

- 0 .4491 741 31 D. 

09 

481 

-0 . 3299 2349 SD 

00 

482 

-0 .432 10 175 7D 

00 

483 

-0 Of? A?0D 

00 

48^ 

0 .1131386000 

01 

485 

0 .611106 35 5 D 

00 

486 

0 ^^ 7^5 ,0(^71 p 

00 

487 

-0 .5755 3869 2D 

00 

49 0 

O.IOOOOOOOOD 

0? 

49 X 

O.IOOOOOOOOD 

0*> 

49 2 

0 .100000000D 

02 

49 3 

0 . 10 00 OOOOOD- 

-0 3 

49 4 

0 . 1OOOO0OOOD-O3 

49 5 

0 . 20 09 09 00 0 D- 

-01 

49 6 

0 .500009900D 

01 

49 7 

0 .POOOOOOOOD 

00 

501 

0 .6 160 00 GOOD 

06 

50 2 

0 .925000000D 

06 

50 3 

0 .56 500 9 900 D 

06 

50 4 

0.4 80000000D 

08 

50 5 

0 .549000000D 

0« 

506 

9 .51^0000000 

08 

507 

0 .860000000D 

08 

508 

0 . 3300 OOOOOD 

03 

599 

O.IOOOOOOOOD 

IT 

510 

0 . 1 .1 0 0 0 0 0 0 0 D 

06 

10 0 .3 

0.2000 OOOOOD 

01 

10 1 0 

0 .100000000D 

01 . 

1011 

0 .200000000D 

01 

1014 

0 .300000000D 

01 

1015 

0 . 30 00 OOOOOD 

01 

1017 

0 .500000000D 

01 

1013 

0 .100000000D 

01 

1019 

O.IOOOOOOOOD 

01 

10 2 3 

0.100000000D 

01 

10 30 

-0 . 100090 000 D 

0 3 

10 31 

0.1. 10000000 D 

0 ? 

10 32 

0 .200000000D 

01 

10 35 

0 . 1000000 00 D 

01 

10.36 

0 .609OOOOO9D 

01 

10 39 

0 .300000000D 

01 

104 0 

0 .2090D0000D 

01 

10 4 1 

O.IOOOOOOOOD 

02 

10 4 2 

0 .150000000D 

02 

104 3 

0 .210090000D 

0? 

19 4 4 

0 . 3QO00O0O0D 

01 

10 4 3 

0 . 10 00 00 000 D 

0.1 

1051 

0.500000000D 

01 

10 5? 

0.200000000D 

01 

1053 

0 .500000000D 

0 .? 

10 54 

0 . 17000000 0D 

0? 

10 55 

O.IOOOOOOOOD 

01. 

1050 

0 . 10 00 00 000 D 

01 

10 6 1 

0 .110000000D 

0 2 

106 2 

0.200000000D 

01 

1063 

0 .200000000D 

01 

10 6 4 

0 .100000000D 

01 

1055 

O.IOOOOOOOOD 

0? 

1066 

o .loooooooon 

01 

106 7 

o.iooooooood 

03 

106 3 

0 .6000 OOOOOD 

01 

1069 

0 . 3000 OOOOOD 

01 

10 70 

0 .950000000D 

02 

1071 

0.600900000D 

01 

10 75 

0 .600000000D 

01 

10 77 

0 . 200000000D 

01 

108.3 

0 .110000000D 

0? 

10 88 

0.120000000D 

02 

10 89 

0.12 3456 70QD 

07 





NO-MIDCOURSE APPROACH CONDITIONS 
RCA INC TCA C3 TRACKING 

(KM) (DEG) (HRS) (K2/S2 ) JTCOPRGS 

40014. 175.77 10 3.33 0.76 82 86.10000? 

MIDCOURSE CORRECTION ANALYSIS 
ANCHOR VECTOR EPOCH 
15. HR 14. MIN 43.26SSEC 
FIXED TIME OF ARRIVAL GniDANCE 


MC TIME MAG. 
JO (HR) (M/SEC) 

1 2.00 45.13 

2 2.00 57.11 

4.00 5 8.14 


FUEL BURN R.A. 
(KG) TIME (DEG) 
(MIN) 

6.72 5.09 14.22 

8.48 6.61 1.58 

8.63 6.74 26.55 


DEC. SPIN- MIN 
(DEG). SUN SUM 

(DEG) (DEG) 
26.89 57.69 4.?f 

16.39 71.83 15.05 

28. IS 46.76 2.2' 

59 


ATT TRACKING 
CHNG VISIBILITY 
(DEG) JTCOFRGS 

148.9 78300000 
162.4 78300000 

137.9 9 3200000 


O 



4 

4.00 

66.85 

9.91 

7.91 

16 .41 

22.0 5 

56,87 

5.53 

14*\9 

93700000 

5 

6.00 

67.64 

10.02 

8.02 

32.69 

2 8.51 

41.42 

1.31 

1 ■»*> . 5 

36000007 

6 

6.00 

75.06 

11.10 

9.07 

2 3.99 

24.06 

49 .70 

3. 55 

140.3 

3600000° 

7 

8.00 

75.6 7 

11.19 

9.16 

35.61 

2 3.61 

33.05 

1.43 

179 . 1 

64000004 

8 

8.00 

82.42 

12.17 

10.16 

28.37 

25.07 

45.74 

2.65 

1 *>6 . 7 

64000004 

9 

10 .00 

82.90 

12 . 2 4 

10.23 

39 .40 

2«.61 

75.67 

1.34 

136! 6 

31000773 

10 

10 .00 

89.23 

13.16 

11.22 

32. 35 

25.65 

42.12 

2.20 

13 7.1 

31000777 


EXEC ERRORS j 0.10 M/S RES. 0.070 PPOP. 0.700 DEG POINT. ( ALL 1 STOMA) 
KISS-VECTOR ERROR STATISTICS 95% PROBABILITY 


MC 

TIME 

TFUSL 

- PERI LOME IMCLIN TFLITE VI N FIN 

VC I PC 

MC2 

FUEL? 


(HR) 

(XG) 

(KM) 

(DEG) 

(SE 

r) (M/S) 

(M/S) 

(M/S) 

(KG) 

1 

2.00 

3.87 

531.1 

16 . 373 

69 4 

.7 3 

.04 

.1 . °° 

15.66 

2.73 

2 

2.00 

4.6 4 

6 34.4 

• 19.7R7 

344 

. 3 3 

.95 

2.77 

13.0? 

2.67 

3 

4.00 

4.09 

550 .5 

16.472 

lid 

.8 2 

.99 

1.9 4 

15 .71 

2.7 7 

4 

4.00 

3.5 8 

6 16 . 8 

18.0 72 

836 

. 3 3 

. 38 

2.30 

16.99 

2.45 

5 

6.00 

4.58 

569 .9 

16 .259 

759 

.0 2 

.99 

2.0 3 

15 . 36 

2.21 

6 

6.00 

4.22 

622.7 

17.457 

376 

.1 3 

.71 

2.45 

16 . 37 

2. 35 

7 

8.00 

5.35 

5 89.6 

16.125 

80 3 

.9 3 

.00 

2.90 

15.6? 

9 . 9 A 

" 8 

3.00 

4.99 

6 35.8 

17. 122 

923 

.3 3 

.76 

2.59 

16 . 40 

2. 35 

9 

10 .00 

6.09 

609 .2 

16 .029 

859 

.4 3 

.03 

? . 

15.88 

2.77 

10 

10 .00 

5.73 

651.6 

16 .904 

973 

.0 3 

.34 

2 . ii 

16.58 

i. I** 




ORBIT 

IN SE PT ION INFO RMATION 




MC 

RETRO 

RETRO 

ATT 

SPIN- 

MIN 

tpik 

T n IM 

TOTAL TRACKING 

MO 

TIKE 

R. A. 

DEC. 

SUM 

SUN 

VEL 

FUEL 

FUEL 

r VTSTBT 

LITY 


(HR) 

(DKG) 

(DKG) 

(DEG) (PEG) ( 

M /S.8C) 

(XG) 

(XG) 

JTCO^RGS 

1 

110.00 

133.78 

-79.57 

75 . 70 

26.24 

29 .46 

3. 

7 p 9 . 

9 4 90000 73 3 

7 

111.00 

138.51 

-29 .52 

74.95 

45 .06 

23.95 

7 

60 n. 

03 00000777 

3 

110.00 

138.76 

-29 .53 

75.16 

20.73 

24.77 

9 # 

6 4 n . 

7R 00000773 

4 

111.00 

13 8.4 5 

-29.46 

74.38 

31.0 9 

30.40 

2 . 

70 I 0 . 

11 0000Q737 

5 

110.00 

138. 73 

-29 .53 

75.14 

18. 75 

20.53 

2. 

22 

74 00000333 

6 

111.00 

138.42 

-29.50 

74.88 

26.04 

1.6 . 7 A 

1 # 

80 p, 

90 00000377 

7 

110.00 

138.69 

-29 .52 

75.11 

16 . 9 3 

17.79 

1 . 

86 17. 

05 00000378 

8 

111.30 

133. 37 

-29.53 

74.36 

2 3.26 

13.45 

1 # 

A A 11. 

61 00000377 

9 

110.00 

13 8.6 4 

-29 .52 

75.03 

15.95 

14.25 

x # 

^9 11 . 

77 00000777 

10 

111.00 

138. 31 

-79 .54 

74.83 

21.44 

10.33 

1. 

71 V. 

76 00000737 


WO MORE DATA, RUM TERMINATED 
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FIGURE 5.2 

MIDCOURSE ANALYSIS SAMPLE CASE 
FIXED ATTITUDE SCAN 

INPUT DATA 

20 6 . 

2 3 15. 

50 S. 

5 3 15. 

40 1051.8313 
43-10.141568 
10 4 4 3. 

1051 2. 

426 .045 
474 5. 


MAESTRO OUTPUT 


DATA CARDS FOR CASE 1 

23 0.150000000D 02 24 0.100000000D 01 25 

50 0 , 60000 OOOOD 01 51 0.100000000D 02 5? 

53 . 0 . 150000000D 02 54 0.140000000D 02 55 

40 0 . 105 13 3130D 04 41 0.593204330D 04 42 

43 -0 . 10 14 156 80 D 02 44 0 . 30 30 210 49 D 01 45 

1044 0.300000000D 01 479 0 . 100000000D-01 1077 

1051 0. 200000 00 0D 01 47 0.520000000D 02 43 

426 0.450000000D-01 1041 0.300000000D 01 1042 

474 0.500000000D 01 0 0.0 0 

LUNAR FIELD t NMOD=*** MMOD=*** 

INITIAL JULIAN DATE 2441844.1352 


INPUT COMMON 


A 

Q.540000000D 

06 

5 

0 . 100000009D 

01 

6 

0.100000000D 

01 

8 

0 . 100000 OOOD- 

-10 

10 

0 .100000000D 

21 

20 

0 . 60 OOOOOOOD 

01 

21 

0.100000000D 

02 

22 

0 . 19 7 3 0 0 0 0 0 D 

04 

23 

0.1500000000 

07 

24 

0.10 0000000 D 

01 

25 

0 . 325000000D 

01 

37 

0 . 41 Q 4419570 

05 

38 

0 . 3 33 390000D 

0 3 

40 

0 . 105183130D 

0 4 

41 

0. 597204330D 

04 

42 

0 .260425590D 

04 

43 

-0.10 14 156 80D 

0? 

A A 

0 . 30 30 77.0 49 D 

01 

45 

-0 .280736700D 

01 

46 

0.4134413520 

05 

47 

0.520000000D 

02 

48 

0 .420000000D 

02 

50 

0 ,6000000000 

01 

51 

0. 100000000D 

02 

52 

0 . 19 7 30000 OD 

04 

53 

0.150000000D 

02 

54 

0.140000000D 

02 

55 

0.4 326 80000D 

02 

100 

0 .216 855 300 D 

01 

101 

0.7243354000 

06 

10 2 

0 . 39 860 3200D 

06 

10 3 

0.4291551500 

05 

10 4 

0 . 126 ' 7 10600D 

09 

105 

0.3791369000 

08 

109 

0 , 132715450D 

12 

110 

0.4907 777 90D 

04 

112 

0.2 329 81600D 

04 

113 

0.60986300D 04 

114. 

0.637316503 D 

04 


0.7250000000 01 
0.197300000D 04 
0.4326000000 02 
0 . 2 6 0 4 .2 5 5 9 0 D 04 
-0.2 80 72.6700 D 01 
-0. 7000000000 01 
0.470000000U 02 
0.300000000D 01 

n.o 


21 10 . 

24 1. 

51 10. 

54 14 . 

41 5932.0433 
44 3.03071049 
479 .01 
47 52. 

1041 3. 


22 19 73. 

25 

52 19 73. 

55 47.260 

42 2604.7559 
45-2.307367 

1077 - 3 . 

43 42. 

104 7 3, 
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115 

0.3409530000 

04 

116 

0.714220000O 

05 

117 

118 

0 .2548400000 

05 

119 

0 .249 83Q0Q0D 

05 

120 

121 

0 . 706000000D 

OS 

122 

0. 17.33000000 

04 

176 

m 

n *7 0 n aa/3A 

.9 * 

1 O O 

A 1 7COC/»OOOn. 

• 

.n -» 

1 O A 

170 

0. 360000000D 

06 

171 

o.ioooooooon 

20 

.130 

181 

3 . i oqonnooon 

04 

19 7 

0 .470000000D- 

-0 4 

19° 

201 

0.900000000D 

03 

202 

0.2000000000 

0 5 

271 

222 

0.200000000D 

02 

260 

0.5337340000 

07 

261 

262 

0.7669920000 

02 

280 

-0.9735000000 

04 

2 81 

321 

0.900000000D 

0 3 

322 

0.2000000000 

05 

331 

332 

0.200000000D 

02 

350 

0.7404000000- 

-01 

351 

357 

0.120200000D- 

■01 

360 

0.353200000D 

01 

361 

380 

0. 100000000D 

21 

381 

0. IOOOOOOOOD 

21 

332 

383 

0. .100000000D 

21 

384 

0 .IOOOOOOOOD 

21 

335 

407 

0,6800000000 

00 

408 

0.1360000000 

09 

409 

410 

0.4835 8326 3D 

00 

411 

0.3255960280 

00 

412 

413 

0 .2599739 4 ID 

0 1 

414 

0.3708619900 

01 

415 

416 

0.494202256D 

0.1 

417 

0.5049821380 

01 

420 

421 

0.116500000D 

0 3 

479 

0 .3960000000 

06 

42 3 

42 6 

0.4 5 0000000 D- 

■01 

433 

0.1356000000 

05 

434 

435 

0 . 700000000D 

00 

4 36 

0.2000000000- 

-01 

437 

438 

0 . 300 OOOOOOD- 

-0 3 

4 39 

0.3600000000 

0 5 

440 

441 

0 .2260000000 

03 

442 

0 .2825000000 

0 3 

44 3 

444 

0.2B3800000D 

04 

445 

0.1314361760 

01 

4 AC, 

44? 

0.600000000D 

04 

443 

0 .6000000000 

04 

450 

460 

0.180000000O 

05 

470 

0 .600000000D 

01 

471 

472 

0.204000000D 

02 

473 

0 . 12 3377000D 

07 

474 

475 

0.200000000D 

02 

478 

0.7200000000 

04 

499 

4 80 

-0.449174131D 

00 

4 31 

-0 . 72992349 3D 

00 

4 37 

483 

-0 .6186844700 

00 

434 

0.113138600D 

01 

485 

486 

0.677 3 6367 ID 

00 

437 

-0 .5755386970 

00 

40Q 

491 

o.ioooooooon 

0 2 

49 2 

0. IOOOOOOOOD 

02 

49 3 

49 4 

0.100000000D- 

-0 3 

49 5 

0.7000000000- 

-01 

496 

49 7 

0.200000000D 

00 

501 

0.6160000000 

06 

502 

50 3 

0.565000000D 

06 

50 4 

0 .4300000000 

08 

505 

506 

0.510000000D 

OR 

50 7 

0. 860000000D 

0 3 

503 

509 

o.ioooooooon 

11 

510 

o.nooooooon 

06 

1003 

1010 

o.ioooooooon 

01 

1011 

0.200000000D 

01 

101 4 

1015 

o.ioooooooon 

01 

1017 

0.5000000000 

01 

10 1 3 

1019 

o.ioooooooon 

01 

1073 

o.ioooooooon 

01 

10 30 

10 31 

o.nooooooon 

02 

10 32 

0 . 200000000D 

01 

10 35 

10 36 

0.600000000D 

01 

10 39 

0.3000000000 

01 

10 40 

10 41 

0. 300000000D 

01 

1042 

0. 300000000D 

01 

1043 

1044 

0. 300000000D 

01 

10 4 8 

o.ioooooooon 

01 

10 51 

10 5 2 

0. 200000000D 

01 

10 5 3 

0.5000000000 

02 

10 54 

1055 

0.100000000D 

01 

10 6 0 

o.ioooooooon 

01 

1061 

1062 

0.200000000D 

01 

10 6 3 

0.2000000000 

01 

1064 

1065 

0. iooooooood 

02 

1066 

o.ioooooooon 

01 

1067 

1068 

0 . 600000000D 

01 

1069 

0. 300000000D 

01 

10 70 

10 71 

0.600000000D 

01 

10 75 

0.6000000000 

01 

1077 

10 83 

o.nooooooon 

02 

10 8 8 

0.1200000000 

07 

10 39 


0.57S050000D 05 
0.^345000000 04 
0 .' 7 29 , ’i:?3SlD-04 

A A/*/- 1 r A A A 

• “* 

O.13OOO0OOOD 05 

0. '>3.00000000 00 

0.1900000000 o? 

0. 70S °9?000D 02 

-0 .9785000000 04 

0.1900000000 O’ 

0.1707000000-01 

0.3537000000 01 

0.3000000000 21 

0.1000000000 21 

0,4000000000 09 

0.1994775100 01 

0.48 36 77 461.0 01 

0 . 28 38000 00 O 04 
0.^200000000 00 
0.7200000000 04 
0 .-'000000000-01 
0.7597000000 06 
0.705797000D 02 
0 . 116 50 0000 D 03 
0.1000000000 21 
0.1673000000 00 
0.5000000000 01 
0. 1000000000-01 
-0.4 373017570 00 
0.633 106 2550 00 
0 . 3.000000000 07 
0. 100000000D-03 
0.5000000000 01 
0.9750000000 06 
0 .5400000000 09 
0.2300000000 08 
0.2000000000 01 
0 . 300000000D 01 
0 . 3.00000000D 01 
-o. iooooooood oi 

0 . 1.00 000000 D 01 
0. 7000000000 01 
0 . 7100 0000 OO 02 
0.200000000D 01 
0.3.700000000 07 
0.1100000000 0? 
0.3,000000000 01 
0.1000000000 03 
0 .9500000000 0? 
-0. 300000000D 01 
0.1734567000 07 


FIXED- ATTITUDE SCAN AT 2.00 HOURS 

DVM RTA DEC RCA INC TFLT ROPA FCP TCF 

35.000 47.000 37.000 3271.1 133.96 3.03.59 3370. 6.32 15.93 
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35.000 

47.000 

42.000 

3396 .7 

136.24 

103.57 

3643. 

7.23 

23.52 

35.000 

47.000 

47.000 

4749 . 3 

139.66 

103.55 

400 2. 

10 .07 

31.31 

35.000 

52.000 

37.000 

3049.5 

132.41 

108.46 

3268. 

7.61 

12.63 

35.000 

52.000 

42.000 

3673.2 

135.03 

108.45 

355 3. 

6.29 

21. 15 

35.000 

53.000 

47.000 

' 4524. Q 

.133.78 

10 3.43 

3915. 

9.43 

29.55 

35.000 

57.000 

37.000 

2954.5 

132.32 

108. 31 

3224. 

8.28 

11.12 

35.000 

5*7.000 

42.000 

3578. 7 

134.99 

10 3. 32 

3514. 

5.82 

20.03 

35.000 

57.000 

47.000 

4431. 8 

133.79 

108. 31 

3830. 

9.1? 

28.78 

DVM 

RTA 

DEC 

RCA 

IMC 

TFl.T 

ROPA 


TCF 

45.000 

47.000 

37.000 

5900.6 

39.62 

110.37 

4164. 

15 . 59 

38.29 

45.000 

47.000 

42.000 

5441.9 

45.42 

110.60 

4034. 

14.6 4 

35.88 

45.000 

47.000 

47.000 

4823.7 

5 3. 79 

110.23 

3847. 

13.0 5 

31.96 

45.000 

52.000 

37.000 

6143.0 

37.77 

110. 30 

4237. 

16.01 

39.48 

45.000 

52.000 

42.000 

5633. 3 

4 3. 46 

110.53 

4096. 

15.03 

36.99 

45.000 

52.000 

47.000 

4957.3 

51.68 

110.17 

389?. 

13.41 

32 . 9 1 

45.000 

57.000 

37.000 

6152.3 

37.07 

110.65 

4?41. 

16.01 

39.52 

45.000 

57.000 

42.000 

56 32.9 

42.71 

110.39 

4096. 

15.01 

36.97 

45.000 

57.000 

47.000 

. 4939.7 

50.9 2 

110.0 3 

3888. 

13.34 

32.80 

DVM 

RTA 

DEC 

RCA 

INC 

TFLT 

ROPA 

FCP 

TCF 

55.000 

47.000 

37.000 

18913.4 

20.6 3 

113. 86 

66 86. 

23.12 

59.43 

55.000 

47.000 

42.000 

17693.0 

23.02 

113.58 

6506. 

2 7.05 

58.88 

55.000 

47.000 

47.000 

15951.4 

26.20 

113.19 

6773. 

22.91 

57.89 

55.000 

52.000 

37.000 

19 331.3 

19 .97 

113. 74. 


2 3.12 

59.65 

55.000 

52.000 

42.000 

13115.5 

on c u 

113.47 

6582. 

23.05 

59.11 

55.000 

52.000 

47.000 

16 321. 3 

25.44 

113.10 

6 299 . 

22.91 

58.15 

55.000 

57.000 

37.000 

19 40 2. 1 

19.6? 

113.55 

6 7P9 . 

23.08 

59.69 

55.000 

57.000 

42.000 

18132.6 

21.97 

113. 30 

6 591. 

23.0 3 

59.13 

55.000 

57.000 

47.000 

16 335.2 

25.09 

112.9 3 

6308. 

22.89 

5 8.17 



FIXED-ATTITUDE 

SCAN AT 

4. DO HOTTRS 



D'/M 

RTA 

DEC 

RCA 

INC 

TFLT 

ROPA 

FCP 


35.000 

52.000 

37.000 

12113.3 

165.09 

10 7.93 

6895. 

15 . 44 

55.5? 

35.000 

52.000 

42.000 

12771.6 

16 3.9 9 

107.85 

7171. 

15 . 4 3 

56.41 

35.000 

52.000 

47.000 

1364 7.3 

16 3.2 4 

107. 76 

7549 . 

15.76 

57.48 

35.000 

57.000 

37.000 

1170 7. 1 

165.0 4 

107.73 

6749 . 

15.3 8 

54.97 

35.000 

57.000 

42.000 

12 391. 1 

163.92 

10 7.72 

7032. 

15.3 8 

55.95 

35.000 

57.000 

47.000 

13295.7 

16 3.17 

107.64 

7418. 

15.34 

57.11 

DVM 

RTA 

DEC 

RCA 

INC 

TFLT 

ROPA 

FCP 

TCF 

45.000 

47.000 

37.000 

5349.9 

146.35 

10 8. A 7 

4188. 

13.24 

36.42 

45.000 

47.000 

42.000 

6178. 3 

145.64 

108.44 

449 7. 

14.6 2 

40.91 

45.000 

47.000 

47.000 

7238. 7 

146.02 

10 8. 39 

4883. 

15 . 80 

45.31 

45.000 

52.000 

37.000 

4770.3 

144.40 

10 8.25 

3970. 

11. 87 

32.52 

45.000 

52.000 

42.000 

5623.5 

14 3.9 7 

10 3.25 

429 7. 

13.71 

38.07 

45.000 

52.000 

47.000 

6710.5 

144.63 

10 8.22 

4700 . 

15.23 

43.31 

45.000 

57.000 

37.000 

4371.5 

143.05 

10 8.0 3 

3817. 

10 . 6 8 

29 . 33 

45.000 

57.000 

42.000 

5237.6 

1.42.32 

108.05 

4158. 

12.9 2 

35.79 

45.000 

57.000 

47.000 

6339 .6 

143.78 

108.04 

4572. 

14.74 

41.75 

DVM 

RTA 

DEC 

RCA 

INC 

TFLT 

ROPA 

FCP 

TCF 

55.000 

47.000 

37.000 

2753.9 

80 . 6 3 

109.20 

30 27. 

11.55 

11.51 

55.000 

47.000 

42.000 

3250.5 

88. 81 

109.12 

3253. 

8.19 

17.43 

55.000 

47.000 

47.000 

3366. 7 

9 8.22 

109.00 

3514. 

11.13 

25.04 

55.000 

52.000 

37.000 

2810.0 

72. «3 

109.0 8 

3051. 

11.10 

11.08 

55.000 

52.000 

42.000 

3204.7 

82 . 10 

10 8.9 8 

32 22 m 

8.46 

16.76 

55.000 

52.000 

47.000 

3713.3 

92. 78 

10 o. 85 

3452. 

10.4 7 

23.35 

55.000 

57.000 

37.000 

2.371.3 

6 7.59 
63 

10 8 . 9 1 

3079. 

10.64 

11.46 



Section 6 

MONTE CARLO ANALYSIS MODE 

This mode is used to determine the probability of satisfying various mission 
constraints. The probability is determined by flying many samples along a desired 
mission flight plan with errors applied at each control point. A .typical mission 
flight plan could consist of the following: 

1. Apply .tracking errors to anchor vector and fly to first 
midcourse time. 

2. Determine first midcourse correction, apply errors to the 
correction, and fly to second midcourse execution time. 

3. Determine second midcourse correction, apply errors to 
the correction, and fly to the Moon, 

4. Determine retro motor firing time and attitude. Apply 
errors to retro and burn retro. 

5. Determine trim fuel to satisfy mission constraints. 

The extent of the mission flight plan is determined through the KMONTE flag 

in location 1052, This flag is defined as follows: 

KMONTE = 3 Two midcourse corrections and retro 

(as described above) 

KMONTE = 2 One midcourse and retro 

KMONTE = 1 Retro only 

As the mission progresses in real time, a point will be reached where the errors 
in the tracking data will reach a minimum. At this time, the next correction or 
retro should be determined for the anchor vector and this correction applied to 
each sample. The subsequent maneuvers are determined in the usual manner. 

The KMONTE flag is set negative to use the program in this mode. • 

The error model is discussed in reference 1. Simply, one sigma errors in the 
point angle and velocity are input for both the retro motor and midcourse motor. 
The midcourse motor also includes a .1 meter/sec . resolution error. 

The user also has the capability of retrieving the state at the beginning of each 
sample. When the KMT OUT flag (location 1074) is set to one, the initial state 
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will be output. The output will be position and velocity vectors in the Earth mean 
equator and equinox of 1950. The random number generator will also be output. 

A Monte Carlo analysis can be restarted at any sample by the input of the random 
number generator in location 1054. 

The inputs for the Trajectory Propagation mode and the Midcourse Analysis are 
also preset for this mode. A guide for the use of this mode is presented in 
Table 6.1. This table presents the preset values of the inputs pertinent to the 
Monte Carlo Analysis Mode. 

A sample of the Monte Carlo Analysis is shown in Figure 6. 1. The first part 
of this figure is the input data. For this run the initial state is loaded as position 
and velocity vectors in location 40-45. The initial time and liftoff epoch are in 
locations 50-55 and 20-25, respectively. A covariance matrix is loaded in locations 
56-91. This matrix represents the injection errors of the launch vehicle. KMONTE 
is set to 3 in location 1052, the sample size set to 3 in location 1053, the auxiliary 
output flag is set to 1 in location 1074, and MODE set to 4 in location 1044. 

The first part of the program output consists of a listing of the input cards followed 
by a listing of the input array including preset values. The Monte Carlo analysis 

follows the listing of the input array. This output consists of descriptions of the 

* 

maneuvers and pre-trim orbit for each sample. A map describing the outputs 
is presented at the beginning of the run. The quantities printed are the fuel, its 
mean and standard deviation for midcourse 1, midcourse 2, trim 1, trim 2, 
trim 3, and rema ini ng fuel. The first two trim maneuvers are associated with 
the Hohmann transfer maneuver to correct radius and eccentricity. The third 
trim is the plane change maneuver to correct inclination. The spin axis-sun angle, 
its mean and standard deviation are also presented for midcourse 1, midcourse 2, 
trim 1, trim 2, trim 3, and the retro. The following lines of output present the 
pre-trim lunar orbit conditions, their mean values and their standard deviations. 
The pre-trim orbit conditions printed are apoapsis radius, periapsis radius, 
eccentricity, argument of perigee, inclination, longitude of the ascending node, 
true anomaly and time of retro fire. 
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A statistical summary of the quantities mentioned above is presented at the 
end of the run. This summary contains the mean value, maximum value 
encountered on any sample, minimum value encountered and standard deviation. 
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Table 6.1 


• Monte Carlo Analysis * ** 

MODE = 4 

INITIALIZED INITIALIZED 


LOCATION 

VALUE 

DESCRIPTION 

LOCATION 

VALUE 

DESCRIPTION 

20-25 

_ 

launch epoch 

475 

20 

retro firing range when 

**30-35 

_ 

initial state 



finding retro conditions 

or 



**1044 

- 

Mode set to 4 

40-45 



**1052 

2 

trajectory profile flag 

**56-76 

- 

covariance matrix 

**1053 

50 

sample size 

**435 

.7 

midcourse pointing error 

**1054 

17 

random number generator 

**436 

.02 

midcourse velocity error 
(percent/ 100 ) 

1074 

0 

initial state output 

**437 

.07 

retro pointing error 

1076 

— 

Element set no. to get 
initial state from GTDS 

**438 

.0003 

retro velocity error 



program 



(percent/ 100 ) 

1080 

0 

retro optimization flag 

439 

36000. 

first midcourse time 

1085 

0 

covariance matrix input 

440 

324000 

second midcourse time 



system 

441 

226. 

midcourse ISP 

1089 

1234567 

satellite ID number 

**442 

282.5 

retro ISP 

1090 

0 

director^ display write 

**443 

71.11 

Mass of retro fuel 




**444 

2838 

desired final orbit 
radius 




**446 

59 

desired final orbit 
inclination 




449 

0 

tolerance band or 
inclination correction 




**472 

20.4 

total midcourse fuel 
available 




473 

13.77 

retro drop weight 





* inputs for a midcourse analysis also set 

** required inputs 
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j tOUTE 

CARLO A’vAI 

vs IS 

s 

£ C?.SE 


IMP IT CARDS 








r 

' * 

21 

10 . 


22 

19 77. 


21 

in. 

24 

40. 


25 

47. 


SO 

r- 

• 

51 

10 . 


5^ 

1977. 


S 2 

in. 

5 4 

49. 


ss 

5 1 . 5 05 


AO 

1020.53175 

41 

59 38. 4123 


42 

260 ° . 207 7 


A 0 - 

l n . 09 0 932 

44 

3.91 °9< 1 41 


AS- 

7.9 709057 


ss 

.21397367 

*> 01 5 7 

.50581777 

n 00 

S°- 

. 505 17035 

n~ 0 1 

S9 

-7 ri q «■» o') *> -> 

n-9 7 <0 

. 1006 °9 °5 

P-01 

G 1 - 

. 710°3 r> 07 

t ^-0 2 

s? 

.1^194174 

n O 1 r *>- 

# ir>//qogQ 

^—0 2 




6 4 

. "■ A 499395 

n-O'J ^5 

.27367741 

P-7 1 

5o 

.'‘4715954 

n ~0 2 

r> ^ 


n i 






£ °- 


r _9 2 cq 

. 76575777 

n — o 2 


. 201 4 0070 

n -01 

71 

. '•'14106 76 

P-0 3 77 

. 75 761 770 

-i _0 *> 


. 04300677 

^^03 

1/ 

. 121171 1 A 

0-0 7 77 

. 3024607? 

D-0 7 




1C. 

. 1.7 'l 9 1829 

n -02 






10 44 

4. 

105 3 

7. 





10 S? 

■3 

' • 

1056 

0. 





10 RO 

1 

a 








MAESTRO OUTPUT 


DATA CARDS FOR CASE 1 


20 

0. 60000 0000D 01 

21 

O.inonOOOOOD 02 

22 

0.1973000000 04 

2 7 

O.150009000D 02 

24 

0.400000000D 02 

25 

0.4 3000000QD 02 

50 

0.6000000000 01 

Si- 

0.1 90000 GOOD 02 

57 

0.197700000D 04 

5 3 

0 . 150000000 D 02 

54 

0.490000000D 02 

55 

0.515850000D 02 

40 

0.102052175D 04 

41 

0. 59 7 04 12 80 D 04 

47 

0.760220730D 04 

4 3 

-0 .100978220D 02 

44 

0.30 1896 1410 01 

45 

-0 .7970.92520D 01 

56 

0.2 179 736 70 D 01 

57 

0 .695813S70D 00 

58 

-0 .505178850D-01 

59 

0.735780220D-02 

60 

0 . 1006 99 350D-0 1 

61 

-0 . 7 10 83 8 830 D-0 2 

€2 

0 . 16 10 4 1740D 01 

63 

-0 . 194489880D-02 

0 

0.0 

64 

0.744000950D-02 

65 

0. 27 16 72 410 D-0 1 

66 

0 .247159540D-02 

67 

0.2496952300 01 

0 

0.0 

0 

0.0 

68 

-0 .6172695700-02 

69 

0 .765852330D-02 

70 

0.7814897000-01 

71 

0.2 74 10 6 760 D-0 3 

72 

0,7536. 137000-03 

73 

-0.343006? 70 D-0 3 

74 

0. 131131140 D-0 2 

75 

0 . 1.02460 2 20 D-0 3 

0 

0.0 

76 

0.1309 18290 D-0 2 

0 

0.0 

0 

0.0 

10 4 4 

0.4000000000 01 

10 53 

0.7000000000 01 

0 

0.0 

105? 

0 . 3000 DO 000 D 01 

1066 

0.0 

0 

0.0 

10 80 

0 . 10 0 0 0 0 0 0 0 D 01 

0 

0.0 

0 

0.0 

LUNAR 

FIELD, NMOD=*** MMOD= 

*** 



INITIAL JULIAN DATE 2441344. 

1596 



INPUT COMMON 





4 

0.5400000000 06 

5 

0.100000000D 01 

6 

0 . 1.0000000GD 01 

3 

0 . louuoooooo-io 

10 

0 . 100000000 D 21 

20 

0.60009000QD 01 

21 

0 . 100 000 000 D 02 

22 

0 . 19 7300000D 04 

23 

0. 150000000D 0? 

24 

u . 4o00uu000i> 02 

25 

G.433000300D 02 

37 

0 .4184415730 05 

38 

0.333390000D 03 

40 

0 , 10 205 2175 D 04 

41 

0.59 38 412 80D 04 

42 

0.260??07 3 o D 04 

43 

-0.100938220b 02 

44 

0.70 1896 141D 01 

45 

”0.29 309 2520 D 01 

46 

0 . 418441596 D 05 

47 

0.163 3 4860 ID 03 

48 

“C .1554613630 02 

50 

0.600000000D 01 

51 

0 .100000000D 02 



52 o . 3.973000000 o* 

55 0.515 850000D 02 

58 -0 .505 1788500-0 1 

fil -0 .710 83 3830 D-0 2 
64 0 . 7 4400 09 50 D-0 2 

67 0 . 24969 52 30D 0.1 

70 0 . 2 8 14 89 30 0 P-0 1 

7 3 -0.34 3006 270 D-0 3 

76 Q.130 9132900-02 

102 0.39 360 3200D 06 

105 0 . 370 ? 96 °00D 0 8 

112 0.2 329 8.16 OOD 04 

115 0 . .34095 300 OD 04 

118 0.254340000D 05 

121 0.7060000000 06 

127 0 . 70 P820000D-04 

170 0.360000000D 06 

181 0.180000000D 04 

201 0.900000000D 03 

222 0 . 200000000D 02 

262 0 .266 p 3 2000D 02 

321 O.90OOOOOOOD 03 

332 0 . 200000000D 02 

352 0 . 120 20000 OD-O 1 

380 0.100000000D 21 

383 0.100000000D 21 

407 0.680000000D 00 

410 0.4 935 8326 3D 00 

413 0 .2599 789 4 ID 0 1 

416 0.49 420 2256 D 01 

421 0 . 1 16 5 0 0 0 0 0 D 0 3 

433 0 . 135600 000 D 05 

J?6 0 .20000 00 00 D-0 1 

439 0.360000000D 05 

44? 0 .2 92 50000 OP 03 

445 0. 131* 35.' 176 D 01 

448 0.60000000 OD 04 

470 0 .600000000D 01 

473 0 . 12 337700 OP 02 

478 0.7200 00000 p 0 4 

481 -0 .32992349 8D 00 

484 0 . 1133 3 9600 D 01 

4 87 -0.57553969 2D 00 

49 2 0 . 100000000 D 02 

495 0. 200000000 D- 01 

501 0 . 6 16000000 D 06 

504 0 , 4 80000000D 08 

507 0.86 0000 GOOD 08 

510 0 . 110000000P 06 

1011 0.200000000D 01 

1017 0.500000000D 01 

1023 0 . 100000000 D 01 

103? 0 . 200000000 D 01 

10 39 0 . 300000000D 01 

1042 0 . 1500 00000 D 02 

10 4 R 0. 100000000D 01 

1053 0 . 3000000 OOD 01 

1060 0 . 100000000 D 01 

1063 0 . 20 OOQOOOOD 01 

1067 0* 100000000 D 03 


t-, 0. 150000000P03 

56 0 .2139 7 36 70 D 01 

59 0. 785R0320D-02 

6? 0.16 104 1740 D 01 

65 0 .2716724 10D-0.1 

68 -0.617269570D-02 

71 0.2 3*106 750 D-0 3 

74 0.131131140D-02 

100 0 .216 55 5300 D 01 

10 3 0 » 429155 150D 05 

109 0. 1 32717450D i? 

113 0.60996. 36 OOD 04 

116 0.714220000D 05 

119 0 .249 830000D 05 

122 0.173900000D 04 

128 0.175 854900 D-03 

171 0 . 1000000 OOD 20 

197 0 . 47000000 OD-O 4 

202 0 .200000 GOOD 05 

260 0 . 5 3 373 4000O 02 

280 -0 . 9*785000000 04 

322 0.200000000D 05 

350 0.240400000D-01 

360 0. 33200000D 03 

381 n.lOOOOOOOOD 21 

384 n.lOOOOOOOOD 21 

408 0 . 136000000P 09 

411 0.925596023D 00 

414 0 . 370 ^6 199 0 D 03 

417 0. 50 49 821 98 D 01 

422 0.396000000D 06 

434 0 . 720000000D 04 

437 0 .'700000000D-0 1 

440 0.2597000000 06 

44 3 0. '705792000D 02 

446 0 . 1 3 6 5 0 0 0 0 0 D 03 

450 0 . 100 0000 OOD 21 

471 0 . 16 7800000 P 00 

474 0 . 500 000000 D 01 

479 0 . 30 000 00 00 D-0 3 

4 32 -0 . 432.10 I757D 00 

485 0. 611106 355D 00 

490 n.lOOOOOOOOD 02 

493 0 . 300 00 000 OD-O 3 

496 0.500900000D 91 

502 0 ,9250000000 06 

505 0 , 5 4OO00 QOOD 0« 

508 0 . 3.30000000D 08 

1003 0.200000000D 01 

10.14 O.3O0OOOOOOD 01 

1018 n.lOOOOOOOOD 01 

10 30 -0 ♦ 100 00 0000 D 0.1 

1035 n.lOOOOOOOOD 01 

1040 0 . 200000000D 01 

1043 0.210000000D 02 

1051 0.100000000D 02 

1054 0.170000000D 02 

1061 0 . 110000000D 02 

1064 0 . 100000000D 01 

1068 0.60^00000D 01 


54 

n. 1 40000000D 

07 

=»7 

0.685813970D 

00 

60 

0. 1.00639 950O-0 


63 

_ 0 .19 448R390D- 

0? 

66 

0.247159540D- 

0? 

69 

0.965952330D- 

02 

79 

0 . *15 361 3 7 00D- 

03 

75 

0 . 1 0?460?20D- 

03 

10 3 

0 . 7748 35400D 

06 

104 

0. 3267106000 

09 

110 

0 .490277 7900 
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10 70 0.950000000D 02 107! 0 . 6000000000 01 1075 O.60OOOOO99D O 1 

1030 0. 100090000D 01 1033 O.IIOOOOOOOD 02 1083 0 , ! 20000000D 02 

10 89 0.12 3456 700D 0 7 

MONTE CARLO ANALYSIS 

ANCHOR VECTOR EPOCH 
15. HP 49. MIN 5!.SP5SEC 

SAMPLE SIZE = 3 

FIRST MIDCOURSE AT 1Q.00HPS 
SECOND MIDCOURSE AT 72.00HRS 
MIDCOURSE GUIDANCE LAW = 2 

1 SIGMA MIDCOURSE POINTING ERROR = 0.70 DEGREES 

1 SIGMA MIDCOURSE VELOCITY ERROR = 2.00 PERCENT 

1 SIGMA RETRO POINTING ERROR = 0.07 DEGREES 

1 S I G’ 'A RETRO VELOCITY ERROR = 0.0 3 PERCENT 

TRACKING COVARIANCE MATRIX 

0.213P' T D 01 0 , 6 P5 RlD 00-0.5051^^-01 0.78578P-92 0 . ! 0009D-01-0 . 710 94D-9? 
0.-R5A1D 00 0.16 10 4 d 0 1-0 . 1<« 44?r-02 0.74400D-02 0 . '>71 ^t , D- 01 0.*>47i fn-0’ 
“9 . 50 51RP-01-0 . 194 49 p-02 0.24970P 01-0 . 6 I 727D-0 7 0.2S585D-02 0 . *>8149 "'-0 * 
0.7B57PD-02 0.744000-02-0.61727^-02 0 .23^! ID-03 0 . ?*! P-0 3-0 . 3439 ■ . ■ 

0 . 10069D-0 1 0.27167D-01 0.26535P-02 0.353Op-03 0.13l!3D-02 0. 1 0 74f!n-0 , 
-0.710 3 40-0 2 0 . 2 4716D-0 2 0. 7 8149D-01-0. 34301D-03 0.!0246D-03 0.13O97n_07 


OUTPUT FORMAT 


MCI 

MC2 

TRIM1 

TRIM2 

TRIM 3 

FUEL 

FUEL 

FUEL 

FUEL 

FUEL 

*HEL 

REGAINING 

(KG) 

(KG) 

(KG) 

(KG) 

(KG) 

(KG) 

MEAN 

MEAN 

MEAN 

MEAN 

mean 

MEAN 


ilC 2 
r UEL 

TRIM], 

FUEL 

TK.U-12 

FUEL 

Trim 3 
FUEL 

ruEL 

REGAINING 
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SIG f *A 

SIGMA 
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SIGMA 
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TRIM2 
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( DEG) 
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MEAN 
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MEAN 
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TRIM1 

TRIM2 
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SPIN- SUN 
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SPIN-SUN 
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SPIN-SUN 

SPIN— SU V 
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SAMPLE 


SAMPLE 


P PE -TRIM I/TNAR ORBIT 


RA 

(KM) 

RP 

(KM) 

EHSEN 

MEAN 

RA 

MEAN 

RP 

MEAN 

ECCEN 

SIGMA 

RA 

SIGMA 

RP 

SIGMA 

ECCEN 

TRUE 
ANOM 
( DEG) 

RETRO 

FIRE 

(HR) 


MEAN 

TRUE 

(DEG) 

MEAN 

RETRO 

(HR) 


SIGMA 

TRUE 

(DEG) 

SIGMA 

RETRO 

(HR) 
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0.013 

2 . 181 
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0 .019 

0.0 

0.0 

0.0 

139. 196 

11.304 

112.142 

138.196 

11.30 4 

112. 142 

0.0 

0.0 

0.0 

3216.3 

2836.6 

0 . 06.2 74 

32 16 . 3 

2 936.6 

0.06274 

0.0 

0.0 

0.0 

-0 .340 

109 . 348 


-0.340 

109.848 


0.0 

0.0 



y.i2i 

0 . 1 SO 

0 . 0 0 o 

5.651 
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0.012 

3.4 70 

0.048 

0 .006 
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135.399 

69.023 

140.263 

73. 351 

90 .032 

2.082 

62.043 

22.059 

306 9. 7 

2838.5 

0,03913 

3143.0 
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0.05094 

73.3 

1.0 

0.01180 

-0.710 

109 .846 


-0.525 

109 .847 


0.185 

0.004 



AOP 

(DEG) 

IN CLIN 
(DEG) 

LAN 

(DEG) 

MEAN 

AOP 

MEAN • 
IN CL IN 

MEAN 

LAN 

SIGMA 

AOP 

SIGMA 

INCLIN 

SIGMA 

LAN 
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13.534 

0 .0 

0 .0 

0.0 

112.132 

119 . 139 
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112. 132 

119 . 1 39 

67.666 

0.0 

0.0 

0.0 

141.270 

116.523 

- 74 . 74 ? 

141.270 

116. 523 

- 74 . 7*7 

9 .0 

0.0 

0.0 


0 or 0 

« w v 0 

0.0 41 

8.183 

3.699 

0.047 

10 . 959 

0.830 

0.006 

2.675 

111.977 

60.738 

67. 434 

112.055 

89.938 

67.545 

0.077 

29.300 

0.121 

141.467 

116.483 

-74. 599 

141. 36 8 

116 .50 3 

-74.66® 

0.099 

0.020 

0.079 
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SAMPLE 


2 




probability 

OF S/TISFYIMU 


SUFFICIEr*™ :: 

NJ-'NUVjEP. FUEL = 

100.00 PERCENT 



S TAT 

ISTICAL S 

rnr *71 d\/“ 



ME A!! 

\* A. V 

mv 

SI^ f A 

MCI REEL 

5.21? 

0 ( V) 

? . 1 9 

7. 09 

.MCI SPIM-S 

10 p # 0 O 


. <7 7 

5./*. 01 

::c 2 fuel 

0.25 

n , An 

n.oR 

0.17 

MC? SPIM-S 

99 .97 

? 5 1 . ? 1 

11 . 70 

61.1? 

m EIM7 FUEL 

0 . 0/1 

0.00 

0.01 

0.0 7 

TRr ' 1 SFIV-S 

0? . 70 

1 1 ? . 1 4 

67.97 

20 . 7< 

T*'!'? ’"'UEL 

2.7- 

/l. ^3 

7. °7 

0 .67 

TRINE spi::-s 

117.97 

in t 1 o 

111.01 

0.17 

TF.I" i t’UHL 

O.or 

0.05 

0 .00 

0.07 
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°0 . 19 

1 19 . 14 

60.60 

77.6/* 
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U. 0 ° 

1 •» . 5 7 

R. i o 

7.21 
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6 7.51 

^7.6-7 

6 7 , A o 

0.11 


P RE— 1-1 

PIN T.rj'lAR 

ORBIT 


RA 

.7146.9 

7216. 7 

7069 . 7 

60.1 

RP 

2 °39 , P 

? Pj A A * A 

2R76 . 6 

7. A 

EOCEM 

0.05.1.7 

0.0627 

0.0191 

0 .0096 

AOP 

14 1. 4 7 

1 4 1. . 6 7 

INI. i7 

0.76 

IN CLIN 

1.16 . 50 

116 .52 

1 1 5 . z* R 

0.07 

LAM 

-74.62 

-74.5/1 

-74. 75 

O .09 

TRUE 

-0.59 

-0.72 

-0.1/1 

0 . 1 R 

PETP.O TINE ■ 

10 9.85 

109. R 5 

103 . 76 

0 .00 


NO MORE r./.T/’ , RUN TERMINATED 
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Section 7 

MIDCOURSE VERIFICATION MODE 

This mode is a complete numerical integration of the midcourse correction. 

The state is propagated from anchor vector epoch to the desired midcourse 
execution time, then the midcourse motor burn is simulated, and the state 
propagated to closest approach to the Moon, This mode also presents an 
extensive tracking station summary and a history of the midcourse burn. 

The velocity away from each of the visible tracking stations is computed 
during the midcourse motor burn in subroutine DOPLER. This data could be 
compared to real-time telemetry data to determine if a non-nominal burn is 
taking place. Any of the trajectory propagation techniques except the averaging 
and multiconic techniques may be used. 

This mode has the capability of obtaining the initial spacecraft state from a 
previous midcourse analysis. This option is specified by setting the KREAD flag 
to the desired midcourse correction number. If this option is not used, the 
initial state, midcourse execution time, burn time, and attitude must be input. 

Table 7.1 presents the inputs to run this mode plus their preset values. 

A sample of this mode is shown in Figure 7.1, The first part of this figure 
presents the input data. This consists of the initial date in locations 50-55 and 
KREAD set to 1. This means that the initial state will be obtained from the first 
midcourse firing time. 

The MAESTRO output consists of three parts. The first part is the input data. 

The input array is listed after the input cards and includes the preset values. 

The second part presents histories of the thrust, weight, and velocity away from the 
visible tracking stations while the midcourse motor is burning. The third part 
of the program output is a summary of the midcourse correction maneuver. 

Pertinent information is presented about the midcourse maneuver, tracking stations, 
and the state at closest approach to the Moon after the maneuver is applied. 
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LOCATION 


Table 7.1 

Midcourse Verification Analysis * 


INITIALIZED 

VALUE 


MODE = 5 

DESCRIPTION 


20-25 

**30-35,1057 
or 40-45 


47, 48 
**50-55 

— 

200-379 

380 

_ 

408 

4.0095D8 

409 

4.0DS 

410-419 

- 

476 

- 

480-489 

- 

**1036 

- 

1045 

1 

1076 

. - 

1087 

- 

1088 

12 

1089 

1234567 

1090 

- 

1092 

— 


launch epoch 
initial state 

spacecraft attitude (used when 1057 £ 0) 
state epoch 

midcourse motor in engine 1 

ignition time when 1057 = 0 

Telemetry carrier frequency # 1 

Telemetry carrier frequency # 2 

observation site longitudes 

midcourse burntime when 1057=0 

observation site latitudes 

trajectory propagator 

doppler analysis flag 

Element set number when initial state 
from GTDS program 

Element set number when attitude input 
from ADP 

unit number of 1087 read 
satellite ID number 
director's write 

Element set number when state transfered 
to GTDS 
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FIGURE 7.1 
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0.1 

106 f 

o . looooonoon 

09 

106 9 

0.600000000D 

01 

1069 

o. iooooooood 

01 

10 70 

0.950000000D 

02 

1071 

0.6000 9 9000D 

01 

10 75 

0. IOOOOOOOOD 

01 

10 77 

0. IOOOOOOOOD 

0.1 

10 83 

0. I.iOOOOOOOD 

01 

10 8 8 

0 . 120 0 00000 D 

02 

10 89 

0 . 17 34 56 700D 

07 
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REAL TIMS THRUSTING INFORMATION 


ANCHOR VECTOR EPOCH 
15. HRS 14. MIMS 41.25 PSEC 
IGNITION TIME = 2. 00 HRS 


TIME 

THRUST 

WEIGHT 

VELOCITY 

AWAY FROM 

TRACKING STATION 

(MIN) ( 

: NEVJT ON) 

(KG) 

JOBURG 

TANANARI 

CARNARVN 

0 . .17 

53.08 

333. 15 

4.04904 

4.09896 

4.07335 

0.23 

52.79 

332.91 

4.04523 

4.08532 

4.07066 

0.50 

52.49 

332.67 

4.04143 

4.08169 

4.06797 

0.67 

52.19 

332 . 44 

4.03765 

4.07808 

4.06530 

0.8.3 

51.90 

332. 20 

4.03388 

4 . 0 7 4 A 7 

4.06263 

1.00 

51.60 

331.9 7 

4.0 3017 

4.07099 

4.35997 

1.17 

51.30 

331.74 

4.07538 

4.06730 

4.05733 

1. 33 

51.01 

331.51 

4.02265 

4.06.773 

4.05467 

1.50 

50.71 

331.2 8 

4.01 89 3 

4.0601 « 

4.75203 

1.67 

50.41 

331.05 

4.01523 

4.05663 

4.04940 

1. S3 

50. 12 

330. 8 ? 

4.01153 

4.05310 

4.04678 

2.00 

49 . 82 

330.60 

4.00 7 85 

4.04959 

4.04 41.6 

2.17 

49 .52 

330 . 38 

4.00418 

4.04607 

4.04155 

2.33 

49.23 

330.16 

4.00053 

4.04258 

4.03395 

2.50 

4 8.93 

329.9 3 

3.996 39 

4. 0 3909 

4.0 36 76 

2.67 

48.63 

329 .71 

3.99826 

4.03562 

4.08377 

2.83 

4 8.34 

329.50 

3.99964 

4.0 7216 

4.07120 

• 

o 

o 

4 8.04 

329.28 

3.9960 7 

4.0787? 

4.07863 

3. 17 

47.74 

329 .06 

7.99244 

4 . 0 7 5 ? R 

4.02607 

3.33 

47.45 

328. 85 

3.9 7896 

4.Q71 86 

4.02351 

3.50 

47. 15 

328.64 

3.9 75 30 

4.01945 

4.0^097 

3.67 

46.85 

328.42 

3.0 ’7174 

4.07505 

4.01843 

3.83 

46.56 

32 8 . 2 1 

3.96920 

4.01166 

4.01590 

4.00 

46.26 

32 8.01 

3.96467 

4.00 029 

4.01339 

4.17 

45.9 6 

327. 90 

3.961 16 

4.00492 

4.01096 

4.33 

45.6 7 

327.59 

3.9 5 755 

4,01.57 

4.00336 

4.50 

45.37 

32 7. 39 

3.95416 

3 . °9 823 

4.00595 

4.67 

45.09 

327.18 

3.95069 

3.99491 

4.00337 

4.83 

44 . 7 8 

' 326 .99 

3.94722 

3.99159 

4.00009 

5. 00 

44 . 4B 

326.78 

3.94377 

3. 933?9 

9.99842 

5.09 

4 4.32 

326.67 

3.94199 

3. 99640 

8.99707 
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MIDcrv’ESE VE RI 7*1 CAT I ON MODE 


ANCHOR VECTOR EPOCH 
6./10. /1.9 73. 

15 . HE \A . MIN 4 3.26 85EO 


MIDCOURSE CORRECTION SUMMARY 

TIMS MACNITFJDE RT.ASC. DEC. FURL BURN TIME 

(HR) (M/SEC) (DEC) (DEC) (K^) (rI vr) 

2.00 45,13 1.4.22 26.39 6.7? ^‘ 09 j 

ATTITUDE ATTITUDE SPIN AXIS 
MANEUVER FUEL SUN ANCLE 
(DEC) (KC) (DEC) 

1.43. P7 14R.RR0 57.69 


TRACKING station summary 


STATION 

LOCAL TIME 

elevation 

AZIMUTH 


(HRS) 

ANGLE 

ANCLE 



(DEC) 

(DEC) 

JOHANNESBURG 

18.0 2 

64.21 

93.94 

TANANARIVE 

20.0 2 

78.09 

155.51 

CARNARVON 

0.02 

25.76 

-106.93 

ORRORAL 

2.0 2 

-0.57 

-116 . 79 

FAT. RBANKS 

7.02 

-56 . 10 

-77.05 

GREENBELT 

11.02 

-52.0/1 

00 . 71 

SANTIAGO 

12 . 0 2 

<r 

• 

cc 

r- ! 

1 

130. 74 


END CONDITIONS AT CLOSEST APPROACH TO MOON 

TIME RADIUS INCLIM C3 

(HR) (KM) (DEG) (KM2/SEC?) 

110.01 3020.57 115.13 0.6159 

NO MORE DATA, RUN TERMINATED 
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Section 8 

ORBIT LIFETIME MODE 


This mode is used to integrate the spacecraft while in orbit about the target 
planet. There are no special analysis techniques associated with this mode. The 
MAESTRO trajectory propagation mode is used with this mode. The orbit lifetime 
was included as a separate mode because many of the trajectory propagation inputs 
that are preset for the trajectory propagation mode must be modified for the orbit 
lifetime mode. The trajectory propagation inputs modified for this mode are. 
shown in Table 8.1 

This mode is preset to fly the lunar orbit using the averaging technique. Thus, 
the compute interval and print table are set to appropriately large values (two 
days). If the user changes the trajectory propagator to a scheme other than the 
averaging, he must also change the compute interval. The multiconic method 
should not be used with this mode. 

A sample of tins mode is presented in Figure 8.1, The inputs used in this run 
were the orbital elements in locations 30-35 and the initial date in locations 50-55. 

The flight time was set to 20 days for tins run. Normally it is 200 days. The 
shadow flag in 1049 and the flag in location 1018 to determine averaged elements 
from input osculating elements are set. The output first consists of listings of 
the input cards array and the entire input with the preset values. The trajectory 
output shown is obtained from the auxiliary output unit. The KOUT9 flag in 
location 1058 is preset to 9; thus, the output is obtained from the output unit nine. 

This output is much more condensed and easier to read for a lifetime case than the 
normal unit 6 output. The unit 6 output is not lost and is available if desired. The 
output in the condensed form is labeled. The only output which might cause confusion 
is FRAME. In this column two integer values are printed. The first integer corres- 
ponds to the planet number while the second integer defines the coordinate system. The 
coordinate system number corresponds to the definition in locations 1039 and 1040. 

The output times shown in unit 9 are shifted from the values in the print table 
because the averaged element start-up procedure was used. 
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Table 8.1 

Lunar Lifetime Mode * 
MODE = 6 


INITIALIZED 


LOCATION 

VALUE 

DESCRIPTION 

4 

. 172 8D8 

Final time 

30-35 



or 

- 

initial state 

40-45 



50-55 

- 

Initial date 

170 

1.D20 

First switching time 

180 

. 172 8D6 

First compute time 

450 

1.D20 

Switching time for printing 

460 

. 1728D6 

Printing interval 

1015 

11 

Launch planet 

1018 

0 

Earth oblateness flag 

1019 

4 

Input coordinate system 

1029 

1 

Lunar oblateness flag 

1030 

1 

Output frequency flag 

1032 

0 

Closest approach test flag 

1035 

1 

Lunar gravity model 

1036 

8 

Trajectory propagator scheme 

1039 

5 

Earth output coordinate system 

1040 

2 

Moon output system 

1044 

- 

Mode flag set to 6 

1058 

9 

Auxiliary output unit number 

1068 

6 

Number of ordinates per interval in Quadrature integration 

1069 

3 

Number of intervals in the Gaussian Quadrature integration 

1098 

0 

Flag to obtain average elements from osculating elements. 


* Inputs to trajectory propagator also set, however 
many locations are modified 
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FI GUPS 9.1 


ORBIT LIFETIME SAMPLE CASH 


INPUT CARDS 


30 

2 8.3 9. 

31 

.004 

3.3 0.0 

33 

95. 

34 

116 . 

35-63. 

10/14 

6 . 

59 

6 . 

' 51 10. 

52 

19 73. 

53 

15. 

54 14. 

55 
100 9 

% 

43. 

1. 

4 

. 177 8 

3 07 10 40 3. 


MAESTRO OUTPUT 








DATA CA 

EPS FOR CASE 

1 







30 

0.28380000013 

04 

31 

0.400000000D-0? 

72 

0.0 


33 

0 .950000000D 

02 

3 A 

0. 116000000H 

0 3 

3^ 

-9 .^OOOQOOOP 

0 2 

1044 

o # *ooooo0(nn 

01 

50 

0 . 600000000D 

01 

51 

Q. 1 000O0090n 

02 

52 

0 . 19 7300000D 

0 4 

5 3 

0 . 1.50000000 D 

m 

5 4 

o # 7 *Aoooooon 

o ? 

55 

0 .430000000D 

0 2 

4 

0. 17280n000n 

07 

1049 

0 . 900000000P 

01 

10 9 9 

0 . 100000000D 

01 

0 

0.0 


0 

o.o 


LUNAR FIELD, NMOD= 

3 ; 

F- -♦ 

o 

u 

11 

3 





INITIAL 

JUL I AM DATE 2441844. 

1352 





INPUT COMMON 








A 

0 . 172300000D 

07 

5 

0 . 100000000 D 

01 

6 

0 . 100000000D 

01 

8 

0 . 10000 0090 D- 

-10 

10 

0 . 10 0 0 0 0 0 0 0 D 

21 

2? 

o . i sooooooo n 

0° 

24 

0 . 140000000D 

03 

25 

0. 430000000D 

07 

30 

0 .28 39000000 

04 

31 

0.4000 oo oo on- 

-0 2 

33 

0 . 050000000 D 

0 2 

34 

0 .1760000003 

03 

35 

-0 .6 30000000D 

02 

37 

0.41 844.1.332D 

05 

27 

0 . 9933900000 

09 

46 ' 

0 .41R441353D 

05 

Al 

0 . 131Q75221D 

03 

48 

0 . 1 4 1 0 0 ° 8 3 7 0 

09 

50 

0 . 600000000D 

01 

51 

o.ioooooooon 

03 

52 

0 . .1973900000 

04 

53 

O. 150000000D 

07 

54 

0. 14-00000000 

02 

5 5 

’ 0 .43000*00000 

09 

10 0 

0 .216 855 300 D 

01 

10 1 

0. 3243354000 

06 

13? 

0 . 39 9 70 00 n 

0^ 

10 3 

0 . 4 39 155 150D 

0 5 

104 

0 , 1 2 6 7 1.0 6 0 0 D 

09 

105 

0 . 9 79 1 °6900n 

08 

109 

0. 132 815 450 D 

12 

110 

0 . 490 ? 77 790 D 

0 4 

1 1 9 

0 9 79° 16000 

04 

113 

0 .609 86 3600 D 

04 

114 

0 * 67^81^50 3D 

04 

1 15 

0 .94095 9900D 

04 

116 

0 , 7 1 4 2 2 0 0 0 0 D 

0 5 

117 

0.575050000n 

!)* 

11 Q 

■ 0.9648400000 

05 

119 

0 . 2 49 9 3000 On 

05 

130 

0.6345000000 

0 4 

7 ? 1 

n nn>r , OOOOOOn 

06 

12 2 

0 . 1 7 3 80 0 0 0 0 D 

04 

126 

0.739 3 1 7 36 i n. 

-0 4 

127 

0.70 88.300000- 

-04 

12 8 

0.175 854900D- 

-0 3 

134 

0.2661.699500- 

■0 5 

173 

0 . IQOOOOOoon 

21 

180 

0 . 172 R0Q000D 

06 

19 7 

0.4700000000- 

-0 4 

19 a 

0 . 91.00000000 

00 

201 

O.900000000D 

03 

202 

0.2000000000 

0 5 

2 ? 1 

0 .■’ 90000000D 

07 

222 

0 . 20000000 on 

02 

260 

0.5 33 7° 40000 

02 

25 7 

o . 9668970000 

02 

2 62 

0 .266 992000D 

03 

2 80 

-0.9785000000 

0 4 

2 

-9 .99Q900000O 

04 

321 

0 .900000000D 

03 

222 

0.2000000000 

0 5 

2 7 1 

0 .1900000000 

0 ? 

332 

0.200000000D 

02 

350 

0.2404000000- 

•01 

751 

0 . 1.908000000- 

•91 

352 

0. 120200000D- 

-01 

360 

0 . 3 5 3 7 0 0 0 0 0 D 

01. 

751 

0.959900000^ 

01 

3 80 

0 . 100000 90 OD 

21 

331 

0. 100000000D 

2 1 

7?? 

0 . 1.000000000 

93 

393 

o.ioooooooon 

21 

384 

0.1000000000 

21 

775 

0.1000000000 

9 i 

407 

0.6900000000 

00 

409 

0.136000000p 

09 

-409 

0.4000000000 

09 

410 

0 .4 3389326 30 

00 

411 

0 . 035 5960790 

00 

A 1? 

0.1994775190 

0 1 

413 

0.2599 7894 ID 

01 

414 

0.970 of 19 9 o n 

01 . 

415 

0 . 4 8 9f 7 7 A f 1 0 

01 

416 

0.49420 22560 

01 

417 

0 . 5040 09 1 o 3 d 

01 

.4 9-0 

0 .9O3900000D 

04 

421 

0 . Ilf 5000000 

03 

422 

0 . 396000000D 

06 

A 2 3 

0.6200000090 

00 

433 

0 .1356000000 

0 5 

434 

0 .730000000D 

04 

*3 5 

0 .700000000D 

00 

4 36 

O.200000000D- 

■01 

437 

o .7000000 oo n- 

•01 

43? 

0 . 800000000D- 

■0 3 
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4 3 n 

0 . 74O0OOOO0P 

05 

440 

0 .?59?00000p 

06 

44 1 

9 .7760099'-tnn 

0 ? 

44 7 

0 .703500000D 

03 

443 

0 . 7*1 5 79 '’O 00 n 

0*» 

444 

0.179 99 090 On 

9 A 

44 n 

O . 1314 36 176 n 

01 

4 46 

0 . 3 1 5 590000D 

07 

447 

9 .^onnooooOn 

0 4 

44-3 

n . ^oooooooon 

0 4 

4 50 

0 .190000000^ 

21 

AC, 0 

0 .1 7->o90990n 

06 

470 

o .4900000000 

01 

4 7 1 

0. 1 6 7«no0OOn 

00 

4 7-1 

9 . in /OOOOOOn 

02 

-4 7.1 

0 . 1? 3 3 7 70 0 ft n 

09 

474 

o.soooooooor 

01 

47^ 

0,7099000991 

07 

47 R 

0 . 7200000000 

04 

47Q 

n # ooooooooon. 

-0 7 

4 00 

—0 . 4401 741 

00 

/!R1 

“0 . 379977/OQn 

no 

422 

-0.497191 •?£; 70 

09 

A 97 

-n # c 1 on 04 47-m 

00 

404 

9 , 11313 363313 

01 

495 

0 .61.1 105 755P 

00 

A PA 

7 . 611 36 96 71 n 

00 

4 P7 

-0.575537-990 

00 

490 

0 . 1009 909 90 n 

9~> 

40 1 

9.19 9000990r> 

07 

49 2 

0 . 100000000P 

09 

49 7 

0 . 1 0 9 0 0 0 0 0 0 r '- 

-0 7 

40 4 

9. 190 0 90990 n- 

-0 7 

49 5 

0.2000000000- 

-01 

496 

0 .ROOOOOOOOn 

01 

40 ^ 

o.ioooooooon 

00 

501 

0.6160000000 

05 

502 

0 . n^OOflOOO r> 

06 

50 9 

9 . 6650 0090 On 

9C 

50 4 

7 . 4900010000 

OP 

505 

0 .54990900011 

9 9 

506 

9 . 6199009000 

07 

507 

O.O^OOOOOOOD 

05 

50 3 

0 . i^nnonono^ 

09 

599 

9. 19090000. 9 n 

1 1 

510 

0. 110000000P 

05 

100.3 

9 ..OOOOOOOOOn 

01 

.10 10 

n.i 00 -0 0000 or 

01 

10 u 

o.poonooooon 

01 

10 1 4 

O.3O0OOO9O9P 

9] 

1015 

9 .119999990^ 

00 

10 17 

0.50 0000000 P 

01 

10 19 

0 . 400990000^ 

01 

19 9 3 

0 . 19 0000000 19 

fp 

10 2 9 

0 . 1.10000000 p 

02 

1:0 70 

0. 10099 Onnon 

01 

10 7 1 

0 . 1 100909900 

97 

10 07 

0 . 10000000 0^ 

01 

10 75 

0 . 1 0 9 9 0 9 9 9 0 0 

01 

10 3C 

9 . 0009090009 

rn 

10 39 

o.sooooooonn 

01 

10 40 

0 . onoaoooooo 

m 

7941 

0,1. 09990 noon 

0~> 

10 A 2. 

0 . 1 5 0 0 0 0 0 0 0 R 

0.7 

10 4 5 

0.2100009990 

0 7 

1 0 4 4 

9 .6999099991 

01 

10 A 0 

•7 . 10000000 on 

m 

10 4 9 

0.2099000O0n 

01 

1.9 51 

0 . 1009009991 

00 

105? 

0 .230000090!? 

01 

10 5 7 

9 . 599000900 0 

97 

19 5 4 

9 , 1 109900999 

97 

1055 

o.ioooooooon 

01 

105 9 

0.00 00099990 

01 

19 6 0 

0 . 1909090001 

ni 

10 5 1 

o.iionooooon 

07 

10 6 7 

0 .2000990900 

0.1 

196 *> 

9. 799 9999999 

0 1 

105/- 

0.10 000 0000 P 

01 

1065 

n . 100 000 9990 

07 

1.066 

0.199 9 900 non 

03 

105 7 

7 . lonooOOOOD 

01 

106 0 

O.69O3OO990D 

01 

1.0 6 0 

0. 700000000P 

01, 

10 70 

0,9509000000 

02 

10 71. 

0. 600000000D 

01 

1075 

0. 6000000009 

01 

10 97 
10 9 p 

O.llOOOOOOOD o? 1093 o.ioooooooon 0° 1009 

o.ioooooooon 01 

PAE-B LH'TAR ORBIT TIME HISTORY 
INITIAL Jf-LI/UT PATE = 3/41344. 1 35 

0.177456 10 on 

07 


TRAJECTORY PPOP ACA-TOi 5 = 
ORDINATES “ 6 INTERVALS = 3 

compute interval = 2. on pays 


tike 

SNA 

KCC5U 

TRUK 

AO 1 * 

TNCTjITI 

LAN 

ERA' IF 

SHAH T 

UPSTART 

WITH 7vVE RACED ELE 

MEETS AT T 

= 6784. 4 SE r. 





2338. : 

l 0.004QP7 -179.53 

5 7.1. 433 

1.23. 762 

-59 .743 

MEAN OF 

1 \50 



0.00 

7 73 3. 3 

0.004037 

-170.635 

9 4.6 75 

115 .997 

-P4.Q71 

11 

7 

0.0 

7. OP 

2238. 3 

0.003342 

9 5 # 2 4 7 

34.. 70 3 

1 15 .9 75 

-90.073 

11 

7 

n.o 

a 9 q$ 

2 73 7. .3 

0.004191 

2, A 17 

75.625 

116.026 - 

n 6 .04 2 

1 I 

7 

0.0 

6 . 0 p 

2773.3 

0 .004162 

-30 .071 

67.043 

116. 199 - 

147.117 

11 

7 

n.o 

0. 9 P 

7^30*0 

0.004909 

- 16 9 . 0 3 2 

Cl .99? 

116.375 - 

167.743 

11 

7 

0.0 

10 .0® 

2733.3 

0.904600 

93.029 

60 .097 

11^ . 792 

1.6 5 . 595 

11 

7 

n.o 

12.0 8 

2773.7 

0.004036 

-3.907 

66,047 

1 1C 9 1 0,0 

1 71 . 453 

1 1 

2 

0.0 

14.0° 

3133.7 

0.904432 

-110 . 700 

76 . 670 

11^.9 00 

117,791 

11 

7 

0.124 

1.6.00 

9 0 7 O 9 

< ' J ' % 1 

0 .004497 

146.153 

05,349 

11" # 975 

97 , 70 4 

1 1 

7 

O.7P0 

3 8.0° 

7 0 7 q # 7 

9.904926 

41 . 471 

Q 4 , 473 

1 15 ,959 

4 1 . 4 1 4 

11 

7 

0 . 759 

20 .00 

7737! 4 

O.OO40S 4 

116 .79 7 

107,943 

11*. 1C 56 

•><c TO 0 

p • 

11 

7 

0.49=: 
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Section 9 


POST-INJECTION TRIM MODE 


This mode is used to determine the trim maneuvers required to adjust the post- 
retro orbit to the desired quantities. The desired orbit is defined as a circular 
orbit of the radius input in location 444 and inclination with respect to the 
equator input in location 446. This mode presents two different trim maneuvers. 
The first is a two-impulse plane change maneuver determined from a method 
developed by F.T. SUN. This method is described in reference 1. The second 
maneuver shown is a two-impulse I-Iohmann transfer maneuver. This maneuver 
corrects only the radius and eccentricity. The inclination is not changed. 

To use this mode, the initial spacecraft state is input in the normal manner, the 
final orbit conditions in locations 444 and 446, and the MODE flag set to seven. 

The input in location 449 is used to control the analysis to some extent. If the 
difference between the inclination of the pre-trim orbit and the desired orbit is 
less than the value loaded in 449, the plane change maneuver will not be determined. 
There are no other special inputs for this mode. 

The output of this mode is shown in Figure 9.1. The first part of this figure 
is the input data cards. The second part is the program output. The analysis 
appears after the input array. This output presents a description of the initial, 
transfer, and final orbits in the target planets equator and prime meridian. The 
second part describes the trim impulses. The attitude of the maneuver is in 
the Earth mean equator and the ecliptic of 1950. The times are shown assuming 
the impulses are performed at the earliest possible times. In actuality, the 
maneuver can be applied at the indicated time plus any integer times the orbit 
period. These times are relative to launch epoch. 

The same output is presented for the Hohmann transfer maneuver. The post- 
trim state is presented in the Earth mean equinox and equator of 1950. 
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The precise post-injection trim maneuver is determined by the targeting portion 
of the post-injection trim analysis. This analysis uses the midcourse logic to 
determine the trim motor's firing time, attitude } and burntime which will result 
in some desired post-trim orbit, The desired end conditions could be obtained 
from the post-injection trim analysis discussed above. The midcourse logic 
uses secant partials to determine the motor’s attitude and burntime at some fixed 
ignition time. However, the analysis automatically scans through a range of 
ignition times so that the characteristics as a function of ignition time can be 
readily accessed. 

To use this mode, the MODE flag is set to seven and the targeting flag in location 
1091 is set to 1. The flags pertaining to the midcourse analysis are not automatically 
set, thus the user should exercise some care when using this mode. Table 9.1 
describes the inputs required to run this analysis. The orbital elements at the 
final time defined by location 4 are the dependent variables used in the iteration 
scheme. If the final time is input zero or less than the required motor burntime, 
the final time is automatically set to the burnout time. A finite burn of the motor 
{location 1071) will initiate a two-step iteration process. The first step determines 
the burntime using the impulsive velocity approximation. The impulsive .velocity 
approximation is used as a first guess in the finite burn step except that the ignition 
time is decreased by half the burn time so that the motor burn is centered about 
the desired ignition time. 

A sample of this mode is shown in Figure 9.2. The inputs indicate that three 
firings are desired beginning at 2735. 8 seconds, spaced at 10 minute intervals. 

The final time is input zero to indicate that the elements at burnout are the 
dependent variables. The desired end conditions are input in locations 420 - 422 
and the tolerances on those elements input in locations 490 - 492. A finite burn 
using Cowell's method is indicated through the flag in location 1071. 
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The MAESTRO output is shown in the second part of the figure. The post- 
injection targeting appears after the input common listing. The output presents 
the motor firing conditions as a function of ignition time. A second block 
presents the post-firing orbital elements. A solution for the second firing time 
is not attempted because the radius at the desired ignition time did not fall 
between the apo-apsis or peri-apsis defined by the desired orbit. 
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Table 9.1 

Post-Injection Trim Targeting 
MODE = 7 


INITIALIZED 


LOCATION 

VALUE 

DESCRIPTION 

4 

540000. 

Initial time. Should reset to zero. 

20-25 

- 

liftoff epoch 

30-35 

40-45 

- 

initial state 

50-55 

- 

State epoch 

410-419 

- 

observation site longitude 

420 

- 

desired semi-major axis 

421 

- 

desired inclination 

422 

- 

desired eccentricity 

434 

7200 

Increment in execution time. 

478 

7200 

First execution time. Should reset. 

490 

10. 

Tolerance on semi-major axis. 

491 

10. 

Tolerance on eccentricity. Should reset to small number. 

492 

10. 

Tolerance on inclination. Should reset smaller. 

1015 

3 

Launch planet. Should reset to 11. 

1019 

1 

Input coordinate system 

1030 

0 

Output frequency flag. Set to -1. 

1044 

- 

Mode set to 7. 

1051 

10 

Number of execution times. 

1063 

1 

Targeting law set to 2. 

1064 

1 . 

Number of times to recompute partials. Should be set to 10. 

1066 

0 

Pre-targeting flag should be set to zero. 

1071 

0 

Trajectory propagator used when burning motor. Do not use 6 

1075 

6 

Propagator used to generate partials. Do not use 6. 

1091 

0 

Targeting flag set to i. 
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Figure 9. 2 
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